NASA/CR — 2004-213199/ VOL2 



Numerical, Analytical, Experimental Study 
of Fluid Dynamic Forces in Seals 

Volume 2 — Description of Gas Seal Codes GCYLT and GFACE 


Wilbur Shapiro 

Mechanical Technology, Inc., Latham, New York 


October 2004 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 

has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/CR — 2004-213199/ VOL2 



Numerical, Analytical, Experimental Study 
of Fluid Dynamic Forces in Seals 

Volume 2 — Description of Gas Seal Codes GCYLT and GFACE 


Wilbur Shapiro 

Mechanical Technology, Inc., Latham, New York 


Prepared under Contract NAS3-25644 


National Aeronautics and 
Space Administration 


Glenn Research Center 


October 2004 



Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http:/ / gltrs.grc.nasa.gov 


FOREWORD 
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1.0 INTRODUCTION 


NASA's advanced engine programs are aimed at progressively higher efficiencies, greater 
reliability, and longer life. Recent studies have indicated that significant engine performance 
advantages can be achieved by employing advanced seals 1 lj \ and dramatic life extensions can also 
be achieved. Advanced seals are not only required to control leakage, but are necessary to 
control lubricant and coolant flow, prevent entrance of contamination, inhibit the mixture of 
incompatible fluids, and assist in the control of rotor response. 

Recognizing the importance and need of advanced seals, NASA, in 1990, embarked on a five- 
year program (Contract NAS3-25644) to provide the U.S. aerospace industry with computer codes 
that would facilitate configuration selection and the design and application of advanced seals. 

The program included four principal activities: 

1. Development of a scientific code called SCISEAL, which is a Computational Fluid 
Dynamics (CFD) code capable of producing full three-dimensional flow field information 
for a variety of cylindrical configurations. The code is used to enhance understanding of 
flow phenomena and mechanisms, to predict performance of complex situations, and to 
furnish accuracy standards for the industrial codes. The SCISEAL code also has the 
unique capability to produce stiffness and damping coefficients that are necessary for 
rotordynamic computations. 

2. Generation of industrial codes for expeditious analysis, design, and optimization of 
turbomachinery seals. The industrial codes consist of a series of separate stand-alone 
codes that were integrated by a Knowledge-Based System (KBS). 

3. Production of a KBS that couples the industrial codes with a user friendly Graphical User 
Interface (GUI) that can in the future be integrated with an expert system to assist in seal 
selection and data interpretation and provide design guidance. 

4. Technology transfer via four multiday workshops at NASA facilities where the results of 
the program were presented and information exchanged among suppliers and users of 
advanced seals. A Peer Panel also met at the workshops to provide guidance and 
suggestions to the program. 

This final report has been divided into separate volumes, as follows: 

Volume 1: Executive Summary and Description of Knowledge-Based System 

Volume 2: Description of Gas Seal Codes GCYLT and GFACE 

Volume 3: Description of Spiral-Groove Codes SPIRALG and SPIRALI 

Volume 4: Description of Incompressible Seal Codes ICYL and IFACE 

Volume 5: Description of Seal Dynamics Code DYSEAL and Labyrinth Seal Code KTK 

Volume 6: Description of Scientific CFD Code SCISEAL. 

This volume describes two industrial codes used to determine performance of fluid film gas seals. 
The codes are GCYLT (gas, cylindrical, turbulent) and GFACE (gas, face), which predict 
performance of a variety of film-riding gas seal configurations. 


^Numbers in brackets indicate references located in Section 9.0. 
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Both codes were written for a PC environment using OS/2 as an operating system. The 
FORTRAN codes, however, are amenable to other systems that use FORTRAN 77, as long as 
memory is sufficient. References 2 and 3 provide the details of code implementation. 

GCYLT is used for analyzing cylindrical gas seals. The turbulent version supplants a previous 
laminar version, GCYL. The laminar results are identical to the previous code, but GCYLT includes 
Couette and Poiseuille turbulence when Reynolds numbers dictate the presence of turbulence. The 
code will automatically determine the presence of turbulence throughout the grid field and make the 
necessary adjustments when computing the pressure distribution. Figure 1 shows solid ring 
configurations and Figure 2 shows typical sectored ring configurations that the program analyzes. 
Program capabilities include the following: 

• Varying geometries, as indicated on Figures 1 and 2 

• Variable or constant grid representation. Maximum grid size is 30 grid points in the 
axial direction and 74 grid points in the circumferential direction. Figure 3 shows a 
typical grid network. The circumferential parameter is 0, and the axial parameter is 
Z. The grid points are identified in the axial direction as I and in the circumferential 
direction as J. The extent of I is 1->M, and the extent of J is 1— >N. 

• Specified boundary pressures or periodic boundary conditions in the circumferential 
direction 

• Axial symmetry option 

• Four degrees of freedom, x and y translations and angular displacements about the x 
and y axes through the seal center 

• Determining load as a function of shaft position or determining shaft position to 
satisfy a given load 

• External pressurization (hydrostatic) of inherently compensated orifices, spot recesses 
or full recesses 

• Choice of English or SI units. 

The output of the program includes: 

• Clearance distribution 

• Pressure distribution 

• Leakage along specified flow paths 

• Load and load angle 

• Righting moments 

• Viscous dissipation 

• Cross-coupled, frequency-dependent, stiffness and damping coefficients 

• Plotting routines for the pressure and clearance distribution 

• Critical mass and frequency. 
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Figure 1. GCYLT Hydrodynamic and Hydrostatic Configurations 
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Figure 2. Floating Ring Concept with Jointed Segmented Rings 
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Figure 3. Unwrapped Seal Surface 
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The GFACE code is used for a variety of seals that can be defined in a polar coordinate reference 
frame. Figure 4 shows typical configurations that the program analyzes. The capabilities of the 
program include the following: 

• Varying geometries as indicated on Figure 4 

• Variable or constant grid representation. Maximum grid size is 51 grid points in the 
radial direction and 361 grid points in the circumferential direction. Figure 5 shows 
a typical grid network. The circumferential variable is 0 and the radial variable is R. 
The grid points are identified in the radial direction as I and in the circumferential 
direction as J. The extent of I is 1— »M, and the extent of J is 1— >N. 

• Specified boundary pressures or circumferential periodic boundaries 

• Three degrees of freedom; axial, z translation and a and p rotations about the x and 
y axes through the rotor origin 

• Determining load as a function of shaft position or determining shaft position to 
satisfy a given load 

• External pressurization (hydrostatic) of inherently compensated orifices, spot recesses 
or full recesses 

• Choice of English or SI units. 

The output of the program includes: 

• Clearance distribution 

• Pressure distribution 

• Leakage along specified flow paths 

• Load capacity 

• Righting moments 

• Viscous dissipation 

• Cross-coupled, frequency dependent, stiffness and damping coefficients 

• Plotting routines for the pressure and clearance distribution 
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Figure 4. GFACE Configurations 
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Figure 5. Polar Grid Mesh System 
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2.0 


THEORETICAL DESCRIPTION AND NUMERICAL 
METHODS FOR CODE GCYLT 


2.1 General Theory 

Reynolds equation for turbulent compressible flow for journal bearings is as follows: 


ill 

>h 3 G 

+ A L 

h 3 G 

r 1 ae ^ 

x dQ) 

<k If 

Z dz) 


= 6pco 


d(ph) 

ae 


+ 12p 


diph) 

dt 


( 2 - 1 ) 


The equation is made dimensionless with the following definitions. (Upper case variables are 
dimensionless). 

Z = z/R, H =h/C o , T =t/t o , P =p/p 0 , 


A = 



12|i/? 2 

P c 2 0 

1 o 0 


, G x and G, = turbulence modifiers 


Substituting the dimensionless variables into the turbulent Reynolds equation produces a 
dimensionless equation. 


d 

de 
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d(PH) 
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d(PH) 
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(2-2) 


For steady-state solutions, the time-dependent term on the right-hand side is eliminated except 
for the computation of spring and damping coefficients. 

In the solution methods subsequently described, the Reynolds equation is not applied directly. 
The Reynolds equation represents the divergence of the mass flow at any grid point. The more 
convenient cell method [4] is to conduct a mass balance directly, and not the divergence of the 
mass flow at each point. 


2.2 Formation of Equations for Determining Pressure Distribution 

Solving for the pressure distribution is accomplished by a method [4] that uses a flow balance 
through a cell volume. The perimeter of the cell extends halfway between the grid point and its 
four neighboring points. A typical cell is shown by the dashed lines on Figure 6. The principal 
grid point is at Row i (length direction) and Column j (circumferential direction). For 
convenience of programming, the grid points are numbered for each cell sequentially from 1 to 
9, with grid point 5 being the principal point. The corners of the cell boundaries are also 
numbered from 1 to 4. 
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Figure 7 shows the flow balance through the cell. There are eight flows across the cell boundaries, 
and there can also be a source (or sink) flow into or out of the cell control volume. The reason eight 
flows are used in lieu of four is that it permits discontinuous clearance boundaries at grid lines (such 
as Rayleigh steps) without taking derivatives across a discontinuous boundary. 

The net flow through a cell can be expressed as: 


Qu 
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034 
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(2-3) 


Qj 2 means the mass flow per unit length across the plus side of cell boundary 1-2, etc. 

The Q's are dimensionless mass flows per unit length, except for Q in which is a dimensionless source 
inlet flow. (Primed values of P are absolute pressures; unprimed values are gage pressures). 

In the 9 direction: 


3 dP 


Q = - PH 3 — G 
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In the length or Z direction: 


Q = -PH 3 



AG 
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(2-5) 


where Q is defined as 


0 = 


12 \iGT a q 

2 r ,3 

PoC 0 


( 2 - 6 ) 


An optional flow can enter the cell from an external source, which can be treated as an inherently 
compensated orifice, or a conventional orifice restriction. Inherent compensation presumes the 
orifice area is the surface area of a cylinder circumvented by the hole size and length equal to 
the clearance under the inlet hole. The conventional orifice area is the area of the hole. The 
conventional orifice generally discharges into a recess that allows the flow velocity to dissipate 
into a region of constant pressure. Two types of recesses are permitted; a spot recess, which is 
treated as a source at one grid point, or a recess of finite length in the axial and circumferential 
directions, which is fed by an inlet orifice. 
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Figure 7. Flow-Balance Across Cell 
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Pressures are taken as the average pressure across the boundary. For example: 

P = P ‘i +P ij * i (2-7) 

12 " 2 


and 


dP | _ PjjH P ij ( 2 - 8 ) 

30 12 AG. 

j 


etc. 

The turbulent G factors are dependent upon the Couette and Poiseuille Reynolds numbers which 
are computed at each grid point 151 . 

The Couette Reynolds number is 


CR(£p 

Re = P H 

\iGT 
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(2-9) 


where the subscript c refers to the cell comer point (e.g., for Q^, P c = P t ). The Poiseuille 
Reynolds number is defined as: 

r: = |vp |//>. 


where R* n = 
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The value of P c is the average of the four surrounding grid points, i.e.. 


p + p + p + p 

p _ i-hV i’ +1 j +1 ij ij+1 


The variation of the G factors with Reynolds numbers are shown on Figures 8 and 9. Once the 
Reynolds numbers have been computed at the cell comer points, the appropriate G coefficients 
are obtained from curve fitting routines. 


NAS A/CR— 2004-2 1 3 1 99/VOL2 


12 



G Factor 



Reynolds Number, hrw p/p 

Figure 8. Couette Turbulence Coefficients 


821602-1 


NAS A/CR— 2004-2 1 3 1 99/VOL2 


13 



G Factor 
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Figure 9. Poiseuille Turbulence Coefficients 
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The final values utilized are: 


G x =A/m[G,(R), G,(0] 


( 2 - 11 ) 


G z = Mb[g,(R), G r («,•)] 


( 2 - 12 ) 


where G x (R e ) and G z (R e ) are Couette, and G p are Poiseuille factors. This procedure selects the 
dominant turbulence factor for each of the eight cell flows. 

By substituting the pressures and pressure derivatives (Equations 2-7 and 2-8) into the mass flow 
balance equations (2-3 and 2-4), an equation is derived that is a function of the nine pressures, P, 
through P 9 , and the clearances taken at the cell comer points, H, — H 4 . Each cell comer point film 
thickness is computed in the clearance routine by appropriate values of Z and 0 and is designated 
as HQ, i = 1,4. For example, HQ is the clearance at the cell comer point 1. 

An optional flow can enter the cell from an external source, which is treated as an inherently 
compensated orifice or the usual hole size orifice restriction. Point sources pose numerical 
instability problems, which are circumvented by applying fine grids surrounding the source 
points. Flow through the orifice is given as: 
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A g = ndoH s C o for Inherent Compensation 


ndo 2 
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for Orifice Compensation 
(spot recess or full recess ) 


(2-15) 
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If 



(2-16) 


where P CR 



(2-17) 


P P i , 

Also , i/ > 1.0, = — — — and P s - P R 


(2-18) 


This condition implies backflow through the orifice. 

The primed values indicate absolute pressure (i.e., f> R = P R +1). 

Thus, a system of numerical equations can be derived as a function of nine pressures. There is 
an equation for every grid point. 

/ (P v P 2 , .... p 9 ) ij = 0 (2-19) 


The system is nonlinear since it is dependent upon multiples of P and its derivatives. 

The solution process starts by assuming a pressure distribution, and using Newton-Raphson 
iteration until the functions f converge to zero within a prespecified truncation error. In equation 
form, the iteration process is: 



( 2 - 20 ) 


where the partial derivatives are explicitly determined, e.g., 


^ =/(/•„/>„ ... />, ... - 

f(P i ->P 2 ' *** Pk-e/2 ’ "• P<)\j 


(2-21) 


The actual convergence is not on f, but on P K (new) - P K (old) , for when the difference vanishes, the 
condition that f = 0 is satisfied. 
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2.3 The Column Method Solution of Newton-Raphson Equations 

The column method 161 is used to solve the new pressures in the set of M x N equations defined 
by Equation 2-20. The advantage of the column matrix method is that its inversions are M x M 
rather than M x N so that its use saves computational time. 

The linearized N-R equations may be written in the form: 

C P +E + D P . = R (2-22) 

j j j j-i j j+i j 

For each value of j, Pj is a vector containing the jth column of new pressures, Rj is the right-hand 
side column vector and Cj, Ej and Dj are in general tri-diagonal matrices. 


Case 1 - Pressure Prescribed at Start and End of Pads 

Equations of form 2-22 are written at all points in the grid corresponding to i = 1, 2, -»M and j = 2, 
3, — >N-1 with boundary column vectors Pj and P N prescribed. Look for a solution in the form: 

P =AP+B (2-23) 

j - 1 j j j 

Where Aj is an M x M matrix and Bj is a vector. Use Equation 2-23 to eliminate P^ appearing 
in Equation 2-22. 


(C + E A .) P + E B + D P . . = R (2-24) 

v j j r j j j j j*i j 

Solve Equation 2-24 for Pj to obtain: 

P = -/ D P , +I.(R. - E B) (2-25) 

j j j j*i j k j j r 

Where I, = (C 3 + E, A,)' 1 (M x M matrix) 

Set j = j+1 in Equation 2-23 to obtain 


P 




B 


j*\ 


Compare coefficients in Equations 2-25 and 2-26. 


A . = -l P , B . . = / (R -E B .) 

j+i j j ’ j*i j v j j r 


(2-26) 


(2-27) 


Set A 2 = 0, B 2 = P t . 

Use Equation 2-27 to compute A 3 , A 4 , --, A N and B,, B 4 - B N . 

Since P N is given and all A } and Bj are computed, we may use Equation 2-23 to compute P N 
P P — P 

r N-2’ r N-3’ ’ r 2- 
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Review of General Procedure for Nonperiodic Boundaries 


1 ) Set A 2 = 0 

b 2 = P, 

2) Compute A j+I , B J+I 

Aj+i = -Ij Dj 

B J+I = Ij (Rj - E s Bj) 
j -> 2 , N-l 

where Ij - (Cj + Ej Aj)' 1 

3) Compute Pj 

p ,i = A i p j + b j 

j -> N, 2 


Case 2 - Column Method for Periodic Boundaries 

Pj , Bj , Rj , Zj are vectors. ft = N -1 

For periodic boundaries, the condition is that P } = P N . At the boundary, j = 1, the general 
equation is: 

c, r, . e, /-, + d, r, = ft, (2-28) 

At column N, the equation becomes 

C„, P,, * E„, P„,_, + O v />, = (2-29) 

To satisfy the boundary conditions, a solution is assumed of the form: 

P .. = A. P. + B. + F P M , (2-30) 

2-1 7 7 7 7 N' 

Aj = 0, B { = 0, F ] = 5 ( Kronecker delta matrix ) (2-31) 

Returning to the general equation: 

C P + E P , + D P , = R (2-32) 

7 2 7 7 ‘1 7 7 + l 7 
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Substituting for Pj_, from Equation 2-30, the following results: 


(C + E .A) P + E B + E F P K , 

' j j r j j j j j n' 


+ D P.. ~R 

j y+i j 


(2-33) 


/ 

j 


= (C. + E AX 

v j j j 


(2-34) 


Then, 


P - -T D P 

J J J 7+1 


+ / (R 

j v j 


E. B ) - / E. F. P N , 


(2-35) 


From Equation 2-30: 


P. =A.. P.. + B . 

j j*i j*i j*i 


Fj * i Tv' 


(2-36) 


Comparing Equations 2-35 and 2-36: 


A . , 

7+1 


^ D 7 B 7+i = W ^ ^ 2 ’ ” N ~ X 


(2-37) 


For P N = P 1} we obtain from Equation 2-30: 


P N , = A 


, P 

A"+l 1 


N ' +1 + P N'*1 ^ 


(2-38) 


After rearranging: 


*V - (5 C<V«. p ! + 


(2-39) 


or 


p - Y 

1 N> l N' 


+ Z N , 


(2-40) 


where 


Y * = (8 - Cv ' + ,) _1 ^ 




(8 


A 1 


(2-41) 
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Substituting Equation 2-40 into 2-30 we obtain: 

- V (F v , F, * Z N ,) + B n , + F„y (F w , F, + Z A „) 

" (^A" ^A" + ^N' yj 1*1 + A n , Z n , + B n , + F n , Z n , 
= F v/ _, F, + Z A „, 


where 


Similarly, 


^V'-l + ^A" » ^JV'-l ^Af' ^A" + ®JV' + ^ N 1 Z N < 


^N'- 2 Av'-i (F^_j Fj + Z^,^) + F^/., + F w /_, (F N , Fj + Z„,) 

= (4at/- 1 ^/-l + F V ,.J F^y) Fj + 4^,^ Z^y j + + F^y^ Z^y 


Therefore, in general: 


— Y P + 7 

1 N / -2 M ^N'-l 


Y . = A. Y + F. F v/ 

y- 1 y y y w 

Z. . = A . Z. + fi. + F Z v y 

y-i y y y y at' 


F , = F. , F. + Z . or P - Y P, + Z 

y-I y-I 1 y-l y y 1 y 


Fj = (5 - F^ 1 z x 


(2-42) 

(2-43) 

(2-44) 

(2-45) 

(2-46) 

(2-47) 

(2-48) 


Review of General Procedure for Joined or Periodic Boundaries 


1) Compute A j+1 , B J+]f F J+] 
A j+1 = -Ij Dj 
B J+1 = I. (Rj -Ej Bj) 
j=l, N-J 


E J+ , = ~h Ej Ej 

A, = 0 

h = (Cj + Ej A/ 1 


Bj = 0 

F, = 8 
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2) Compute Y N , Z N 

Y' = (S -Fj\ 

A = (5 -F N y'B N 

3 ) Compute 

Y H = A i Y J + F J K 

j =N-> 2 

Zj-i - Aj Zj + Bj + Fj Z N 

4) Compute P , = (5 - Yj)' 1 Z } 

5) Compute P j = Y } P l + Z j 
j = 2-^ti 

The coefficient matrices Cj, Ej and D J? and the right-hand side vector Rj, are easily formulated. 
Cj contains all the coefficients multiplied by Pj. By examining Equation 2-20, it is seen that for 
any row i and column j that values of C are: 


1J 3p 8 

% 


dP 5 


(2-49) 


C 


M+1J 


3/ s 

3P 2 


Similarly, the coefficient matrix Ej contains the elements: 


and 


E 


3/ s 

dP „ 




3/ E 3/ 5 

a p/ ap 7 


(2-50) 


Z) ... = 


3/s 

ap. 


D 


m+u 


ap. 


D 




3/s 

ap n 


(2-51) 


Rj contains all elements not multiplied by the pressure 


R = ~f, 

j 


(old) 


E 


A 

•A 


p («/</) 

' K 


(2-52) 
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2.4 Film Thickness Distribution with Eccentricity 
and Misalignment (see Figure 10) 

In vector format, the clearance due to eccentricity and misalignment at any angle 0 and at 
distance z from the mid-plane is: 


h = |c o e r - e x i - ej - ai x z'k - $j x z'k] .e r 


(2-53) 


h = C g - e x cos0 - e^sinO + a z‘ sin0 - P^cosO 
= C - (e + pz') cos0 - (e -a z f ) sin0 


(2-54) 


Using dimensionless variables, Equation (2-54) becomes: 

(Z -L/2)R 


H =1 - 


e x 


e - a 


(Z -L/2)R 


COS0 


sin© 


which is set equal to 


(2-55) 


H = 1- (e x + ej cos0 - (e + e a ) sin0 


(2-56) 


where 


0 _ n (Z-LI2)R 
p r 

o 

(Z -LI2)R 

e = a 


(2-57) 


Preloaded Configurations 

Preloaded configurations (see Figure 11) can be modeled by adding an additional eccentricity in 
the x and y directions. 




cos0 

p 

sin0 

p 


(2-58) 
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Keyword 

Variable 

Description 

START 

THETAS 

Pad Start Angle 

PAD ANGLE 

THETAP 

Pad Angle 

PIVOT 

PVT 

Pivot Angle 

PRELOAD 

EPR 

Offset/Clearance 


Figure 11. Preloaded Configuration 
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where 


e p £ = x eccentricity due to preload 
£ p p = y eccentricity due to preload 
0 p = preload angle. 


Rayleigh Step 

The grid network for the Rayleigh step is shown on Figure 12. The boundaries of the step are defined 
by the lower left and upper right corners of the depressed region. Interior grid points include the step 
height in the clearance distribution. 


Axial Taper 

An axial taper is indicated as Figure 13. If Z > Z t then 


H = H + 8 (Z -Z) 


(2-59) 


2.5 Power Loss (Torque) 

Power loss is obtained by integrating the viscous shear forces across the film. From Figure 14, 
a force balance on an element produces: 


dp _ dx 
dz 


but 


X 



Therefore, 


dp p du 2 

dx Gx dz 2 


(2-60) 


(2-61) 


(2-62) 


where G t is the turbulence shear modifier. See Figure 8. 
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Figure 12. Rayleigh-Step 
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r + dr/dz 



d p/S x 


95TR34-V2 


Figure 14. Viscous Power Loss 


Integrating, 


= 2l ^IL Z + C (2-63) 

dz ji dx 1 

U = —L ^E- - + C.z + C 2 (2-64) 

p dx 2 1 2 


The boundary conditions are: 


U = 0 z = 0 c 2 = 0 
U - U when z - h 


Substituting: 


_ dp h 2 
\i dx 2 


C { h 


(2-65) 


( 2 - 66 ) 


NAS A/CR— 2004-2 1 3 1 99/VOL2 


28 




Therefore, 


and 


c = — - !L 1 

1 h [i 2 


(2-67) 


v = G dp 
fl dx 


u 

+ — z 
h 


(2-68) 


dll 1 dp 

dz fi dx 


h 

z ~ — 

2 


U 

h 


(2-69) 


|U du _ dp 
G dz Ox 


h 

IZ ~2 


V_P i 
/i G 


(2-70) 


, , dp h ju U 

t (at z=h) = + — 

dx 2 G h 


(2-71) 


= friction force = jjldA 


II 


^P_ + M_ j?®] 

foo i g h 


n Po C o H dP + p Fco 

m on r nr 


F 2 dQdZ 


,, r> h dP u co F 3 
C R +— 

00 2 ae g ch 


dQdZ 


(2-72) 


FF = 


JJ 


H SP + A £>. r _L 

0 0 2 ae g 6 // 


</0c/Z 


FF 


JJ 


H dP A 

T ae G6H 


dddz 


(2-73) 
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TF 


ff[- 

dP 

A 

+ 

JJ[ 2 

30 

G6H 


dQdZ 


(2-74) 


G, the turbulence modifier for friction, is the minimum of the Couette and Poiseuille values. The 
Couette and Poiseuille Reynolds numbers are computed as before (see Section 2-2). For the 
Couette number, G t , as shown on Figure 8, is obtained. For the Poiseuille flow, G p , as shown 
on Figure 9, is obtained. The value of G is the minimum of G T and G p . 


2.6 Computation of Flows 

The program computes the flow across specified axial and circumferential grid lines. A total of 
four grid lines can be prespecified. 


Circumferential Flow Line (see Figure 15) 

There are three types of points to consider. A point on a grid boundary J = 1 or J = N, and an 
interior point. Also, a flow line on I = M requires special consideration. For each point, a flow 
balance is accomplished through the cell surrounding the point as depicted on Figure 15. 
Consider an interior grid point on an interior grid line (I ^ M). 


Q c (L J ) - Qu + fil4 + Q 12 Qm 


(2-75) 


where 


fin 



(2-76) 



- r„,, 


-p )/ ae. 

‘JJ J 


(2-77) 


p n = n, 


‘V-1 


)/2.0 


(2-78) 


The remaining flow components are similarly computed and Q c (I, J) determined. 
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At J = 1, 3P/30 1 34 is computed by forward difference and is equal to 3P/30 

The pressure P 34 + = P 12 +. 

The clearance H 4 is not a regular grid point clearance and thus is not included in the grid 
clearance array. H 4 is computed as the average of and H i+ij . 

The grid line mass flows are accumulated to obtain the total flow across the grid line. 

Similar procedures are used for computing flows across axial lines (see Figure 16). 

2.7 Frequency-Dependent Spring and Damping Coefficients 

Discretization has been carried out with the use of the cell method 141 , which involves a flow 
balance about each grid point. 


jV QdA = jQ ndS = -JLj(l +P) HdA 


(2-79) 


where Q = the mass flow vector per unit length. 

The equality of the first two terms comes from the divergence theorem. 
In numerical format, the right-hand side becomes: 

-—1(1+/*. ) H A 1 
07 r */ ‘j ‘j J 


where, 


A.. - — (A0. +A0. ,) (AZ.+AZ. ,) 

ij 4 J J ~ l 1 l ~ l/ 


(2-81) 


Generally, a small perturbation analysis is used for determining frequency-dependent spring and 
damping coefficients and solving the complete equation (2-79). A small perturbation analysis, 
however, is generally limited to concentric operation and produces complex expressions for the 
perturbation coefficients. Identical results can be achieved by direct numerical perturbation of 
the difference equations used in the column matrix solution approach. This method, which is 
described below, avoids algebraic error in determining the perturbation coefficients and may be used 
in complex situations where analytical determination of the perturbation coefficients is not feasible. 
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After desired convergence of the Newton-Raphson process has been achieved under steady 
(unperturbed) conditions, the resulting steady-state pressure vectors are denoted as {£} and the 
coefficient matrices as [£'], etc. and, as before, the steady state becomes: 

[c'lp,} + [£% ,} + H{p„} = {«'} < 2 - 82 ) 


The eccentricity components can be perturbed individually by an amount, r|, and the matrix [£)] 
recalculated at the new film thickness (but old pressure distribution, £); then subtract [6>] at the 
old film thickness and divide the difference by r\ to numerically obtain the partial derivative of 
[£ j ] with respect to the eccentricity perturbation. This partial derivative will be denoted by [£>*]. 

Thus, 


[£/*j = [^1^ (2-83) 

n 

The matrices [£ j k ], [6 i,k ] and [k j k } are obtained in a similar manner. Equation (2-78) may now 
formally be differentiated with respect to e k to obtain the expression: 

IP'AhHP'A} = <2 - 84) 


where {f > j ,k } = 3{h J }/3£ k is the zero-frequency stiffness pressure. This expression does not yet 
contain the time-dependent terms found on the right-hand side of Equation (2-75). It is assumed 
that a sinusoidal disturbance is applied to the shaft, such that the clearance and pressure 
derivatives are affected as follows: 


H = e ia ‘; P* = pfe iai (2-85) 

To complete the process, the right-hand side of Equation (2-75) is differentiated with respect to 
£ k and the results added to the right-hand side of Equation (2-80) with d/dt replaced by io. The 
terms to be added to the right-hand side of Equation (2-80) in this manner are -ia [C j ](£j ,k )- 
ia{R J,k } where [C J ] are diagonal matrices whose components are 

C j = H.A. t 2 ' 86 ) 

a r j ij 
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Because a cell can have clearance discontinuities, such as a step, it is advantageous to partition 
the cell into four components as indicated on Figure 7, and then sum the components to obtain 
[C J ]. Thus, Equation (2-82) becomes: 


cj = HC ] A X + HC 2 A 1 + //C/F + //C/l 


(2-87) 


where HC, is the clearance at the corner point 1 of the cell and 


(A0)(AZ) 
A i = — h 


(2-88) 


(A0.)(AZ i . 1 ) 

A. = ; etc. 

* A 


and {R j,k } are column vectors whose components are 

R J y (1+P..) 

' y de k " 


(2-89) 


By combining terms, the final set of linear difference equations for the complex stiffness pressure 
derivatives (Pj k ) are obtained 

[c -'Ip/} +[£ j ]{p/,} l} -[o'lpj.i} (2 - 90) 


where, 


[C*j] = [cJ] * io[c J ] ; = {r w }-,-o{r'’*} ( 2 - 91 > 

The system of equations given by Equation (2-86) is solved by the column method in a directly 
analogous manner to that used in solving the steady-state equation. The principal difference is 
that all the matrix operations are performed using complex arithmetic. Integration of the real part 
of the pressure derivatives yields the stiffness; and the complex parts, when integrated and 
divided by a, yield damping. 


2.8 Critical Mass and Frequency 

A critical mass routine was also added to the turbulent version of GCYLT. 

For incompressible theory, a closed-form solution exists for the condition of neutral stability for 
a two-degree-of-freedom, point mass supported by cross-coupled springs and dampers [7] . The so- 
called critical mass and frequency are obtained, which provides a measure of the stability margin. 
If the mass attributable to a seal or bearing exceeds the critical mass, then an instability may 
occur at the orbital frequency calculated. A similar analysis, including the effects of 
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compressibility, is complicated by the frequency dependence of the coefficients. A solution exists 
when the computed frequency of the point mass is equal to the excitation frequency used to 
compute the stiffness and damping coefficients. The critical mass and the orbital frequency are 
given by the following equations: 


M = 


BED 

E 2 +AED+CD 


(2-92) 


0) 


( AED+E 2 +CD 2 ) 


ED 


where A -K +K , B -D D -D D 

yy xx 7 yx xy xx yy 


(2-93) 


C =K K -K K ; D =D +D 

xx yy xy yx 7 yy xx 


E=D K +D K -D K -D K 

xy yx yx xy xx yy yy xx 


A Newton-Raphson algorithm was utilized. For convergence, the frequency assumed in 
computing stiffness and damping, co a , should equal the frequency computed by the critical mass 
equations, co c . Initially there will be an error, 8. 

co - co = 5 (2-94) 

a c 

To compute the incremental change in co a , the following equations apply: 

6 oW + i5_Aco = 0 (2-95) 

0(0 

a 


but from Equation (2-94) 



3co 3co 

a a 


2 2 


-E E 

'=1 j=i 


0co dK.. 0co dD.. 

u + f_ i j 

dK.. 0(0 dD.. 0co 

y “ y ° 


(2-96) 


The partial derivatives, do) c and 3 cd c . can be computed directly from Equation (2-99). 

dc^ dcOjj 

The derivatives dK S j and c)D IJ are explicitly determined by computing the gas seal stiffnesses at 

dco 2 3 co 2 

two incremental frequencies and computing AK ;j and ATE, respectively; 

Aco 2 Ao ) 2 
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and from Equation (2-95) 


-5°" 

~db~ 


Aco 

a 


d(0 


( 2 - 97 ) 


The new value of 0t) 2 for the next iteration is obtained from: 

CO r - «£" + Aco a ( 2 - 98 ) 

At times, the computed frequency, co c , was insensitive to variations in the assumed frequency, co a . 
In such cases, successive substitution was found to converge very rapidly. 

The critical mass as defined above applies to a two-degree-of-freedom system. It does not 
include cross-coupling moments. Thus, it would not have application to seals with pressure 
gradients where moment cross-coupling may be significant. 


2.9 Force Balance 

Newton-Raphson algorithms are used to produce journal displacements to balance given loads, 
if this option of the code is applied. The requirements are that: 

F X « ~ Fxc = 0 < 2 -"> 

F yr - Fy C = 0 ( 2 - 100 ) 


where, 


Fxr 

= required load in X direction 

Fyr 

= required load in Y direction 

Fxc 

= calculated load in X direction 

p 

17 YC 

= calculated load in Y direction 

but, 



F = F + A F 

1 XC 1 XoU X 



AX + 


dF x 

a y 


AY 


(2-101) 
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( 2 - 102 ) 


F = F + AF - F + 

YC Void y Yold 


dFy 


AX + 


dF x 

dY 


AY 


The partial derivatives are stiffnesses which are obtained by incrementing X and Y displacements. 
Equations (2-101) and (2-102) are substituted into Equations (2-99) and (2-100). The system 
of equations are solved for AX and AY which provides the new displacements for the next 
iteration. The process continues until the calculated load equals the required load within the 
convergence criteria. 
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3.0 SAMPLE PROBLEMS FOR CODE GCYLT 


3.1 Sample Problem 1 - Rayleigh-Step Seal 

A four-pad Rayleigh-step seal, without pressure differential, was analyzed (refer to Figure 17). 
The geometry, operating conditions, and input features were as follows: 

• Number of pads = 4 

• The seal position was prespecified 

• Seal length = 3.852 in.; the symmetry option was used 

• Variable grid was used in the axial and circumferential directions. Since symmetry has 
been applied in the axial direction, the variable grid length equals half the actual length, 
and is equal to 1 .926 in. 

• The grid was made finer at the step boundaries where sharp pressure gradients are 
expected to occur 

• Seal diameter = 1 .9685 in. 

• Step height = 0.00165 in. deep, located at the leading edge of the pad, 5° from the x axis; 
the lower left comer of the step is 0.655 in. from the left axial boundary. The end of the 
step is 70.6° from the x axis, and since symmetry has been invoked, the axial end of the 
step as represented on the grid is 1.926 in. from the inlet end, or at the end of the grid. 

• Specific heat ratio of the lubricant =1.4 

• Gas constant = 250,000 inV(s 2 -°R) 

• Absolute temperature = 530°R 

• Absolute viscosity = 3.0 x 10' 9 lb-s/in 2 

• Eccentricity ratio = 0.4 

• Eccentricity angle = 270° 

• Shaft speed = 70,000 rpm 

• Reference pressure = 200 psig 

• Boundary pressures are all 0 psig, or 200 psia. 

Summary results are shown on Table 1. Figures 18 and 19 show the clearance distribution and 
the pressure distribution produced by the plotting programs. These plots clearly show the highly 
loaded pad, which is pad number 3 (highest pressure level and lowest film thickness level). 
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861599 

Figure 17. Rayleigh- Step Geometry 
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Table 1. Sample Case 1: Raleigh Step Seal 


Journal & Load Position 
Eccentricity 
Eccentricity Angle 
Minimum Film 
Load 

Load Angle 
Power Loss 
Leakage at I = 1 


= 0.4000 

= -90.00° 

= 0.0006015 in. 

= 131.5 lb 

= 58.49° 

= 2.308 hp 

= 0.10191 E-02 lb/s 
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Figure 18. Rayleigh-Step Seal Clearance Distribution 
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Figure 19. Rayleigh-Step Seal Pressure Distribution 
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3.2 Sample Problem 2 - Nongrooved Lobe Seal 

The nongrooved lobe seal is characterized by offset lobes that are joined at their appexes in a 
continuous fashion as opposed to a lobe seal where the lobes are separated by axial grooves. 
Such a seal is depicted on Figure 20; it would be manufactured by a broaching process. For this 
example, a lobe hydrodynamic geometry was combined with external pressurization through 
source points at the mid-plane of the seal. The geometry and operating conditions are as follows: 

• Seal Diameter = 2.25 in. 

• Seal Length = 1.625 in. 

• Seal reference clearance (the clearance prior to preload) = 0.0005 in. 

• Preload on each lobe is 0.5, which means the lobe is eccentric toward the shaft by a 
distance of one-half of the reference clearance 

• Viscosity of the gas = 3 x 10~ 9 lb-s/in. 2 

• Gas constant = 2.5 x 10 5 in 2 /(s 2 -°R) 

• Ambient temperature = 510°R 

• Total number of orifices = 27, 9 in each sector, located at the mid-plane of each sector; 
1 orifice is located at each interior grid point at the mid-plane of the bearing 

• Orifice diameter = 0.015 in. 

• Coefficient of discharge of each orifice = 0.9 

• Supply pressure to the source orifices =120 psig 

• Operating speed = 70,000 rpm 

• Reference pressure = 14.7 psig 

• Pressure along the boundaries = 0 psig. 

An output summary is shown on Table 2. and the clearance and pressure plots are shown on 
Figures 21 and 22. Notice the discontinuities in the clearance distribution because of the lobed 
geometry. The proximity of the source points to each other makes the pressure distribution 
appear as a line source. 
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Table 2. Sample Case 2: Sectored Seal 


Journal & Load Position 
Eccentricity 
Eccentricity Angle 
Minimum Film 
Load 

Load Angle 
Power Loss 
Leakage at I = 1 
Leakage at I = 1 


= 0.00000 
= 0 . 00 ° 

= 0.0002500 in. 

= 0 lb 

= -54.45° 

= 1.251 hp 

= 0.14815 E-03 lb/s 

= 0.14815 E-03 lb/s 
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Figure 21. Clearance Distribution, Sectored Lobe Seal 
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Figure 22. Pressure Distribution , Sectored Lobe Seal 
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3.3 Sample Problem 3 - Three-Lobe Seal 

This problem deals with the hydrodynamic portion of a three-lobe seal where the lobes are separated 
by axial grooves. Figure 23 shows the general geometry and key parameters. The principal 
parameters are the preload and pivot angle. The following are geometry and operating conditions: 

• The position of the seal to satisfy a given load was to be determined 

• International units were invoked 

• Stiffness and damping was to be calculated in two degrees of freedom, x and y, at an 
imposed frequency equal to running speed of 50,000 rpm 

• Number of pads = 3 

• Start of the first pad is at 100°; the pad extent is 100° 

• Pad preload is 50% of the reference clearance, and the preload for the first pad occurs 
150° from the x-axis, which means the preload is in the center of the pad 

• Shaft diameter = 0.0508 m 

• Hydrodynamic length = 0.0254 m 

• Reference clearance or machined pad clearance = 1.27 x 10' 5 m 

• Lubricant viscosity = 2.07 x 10" 5 N-s/m 2 

• Absolute temperature = 283°K 

• Ratio of specific heat of the gas = 1.4 

• Gas constant = 290.32 m 3 /(s 2 -°R) 

• Symmetry is applied in the axial direction 

• Load to be supported = 200.16 N 

• Angle at which the load is applied was 270° from the x-axis 

• Initial eccentricity guess is 0.2; initial displacement angle guess is 270° from the x-axis 

• Shaft speed = 50,000 rpm 

• Reference pressure = 8.274 x 10 5 Pa. 

• Boundary pressures are all 0 gage. 

Summary results including stiffness and damping coefficients are shown on Table 3. Graphical 
representation of the clearance distribution and pressure distribution are shown on Figures 24 and 
25, respectively. Figure 26 shows the pressure distribution viewing along an axial direction. The 
negative pressures are induced by divergent clearance regions. The final eccentricity ratio to 
balance the applied load was 0.22, and the eccentricity angle was 129.42°. 
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Keyword Variable Description 


START 

PADANGLE 

PIVOT 

PRELOAD 


THETAS Pad Start Angle 

THETAP Pad Angle 

PVT Pivot Angle 

EPR Offset/Clearance 


Figure 23. Three-Lobe Seal 
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Table 3. Sample Case 3: 

Three-Lobe Gas Seal 

-JOURNAL & LOAD 

POSITION 



ECCENTRICITY 


= 

0.22087 


ECCENTRICITY , 

ANGLE 


129.42 

DEG 

MINIMUM FILM 


= 

0.0000037 M 

LOAD 


= 

200.2 

N 

LOAD ANGLE 


= 

-90.00 

DEG 

POWER LOSS 


= 

186.2 

W 

LEAKAGE AT I = 

= 1 

=- 

■0. 12904E-04 

KG/S 

-STIFFNESS COEFFICIENTS 




PRINCIPAL X 

KXX 

= 

0. 1459E+09 

N/M 

CROSS-COUPLED 

KXY 

= 

-0.2373E+08 

N/M 

CROSS-COUPLED 

KXA 

= 

0 . 0000E+00 

N/RAD 

CROSS-COUPLED 

KXB 

= 

0 . 000QE+00 

N/RAD 

CROSS-COUPLED 

KYX 

= 

-0.3890E+08 

N/M 

PRINCIPAL Y 

KYY 

= 

0. 1003E+09 

N/M 

CROSS-COUPLED 

KYA 

= 

0 . OOOOE+OO 

N/RAD 

CROSS-COUPLED 

KY8 

= 

O.OOOOE+OO 

N/RAD 

CROSS-COUPLED 

KAX 

= 

0 . OOOOE+OO 

N-M/M 

CROSS-COUPLED 

KAY 

= 

O.OOOOE+OO 

N-M/M 

PRINCIPAL A 

KAA 

= 

0.0000E+O0 

N-M/RAD 

CROSS-COUPLED 

KAB 

= 

0.0000E+00 

N-M/RAD 

CROSS-COUPLED 

KBX 


O.OOOOE+OO 

N-M/M 

CROSS-COUPLED 

KBY 

= 

0.0000E+00 

N-M/M 

CROSS-COUPLED 

KBA 

- 

O.OOOOE+OO 

N-M/RAD 

PRINCIPAL B 

KBB 

- 

O.OOOOE+OO 

N-M/RAD 

-DAMPING COEFFICIENTS 




PRINCIPAL X 

DXX 

= 

9977. 

N-S/M 

CROSS-COUPLED 

DXY 

= 

-0 . 1090E+05 

N-S/M 

CROSS-COUPLED 

DXA 

- 

O.OOOOE+OO 

N-S/M 

CROSS-COUPLED 

DX8 

= 

O.OOOOE+OO 

N-S/RAD 

CROSS-COUPLED 

DYX 

= 

5162. 

N-S/M 

PRINCIPAL Y 

DYY 

= 

0. 1392E+05 

N-S/M 

CROSS-COUPLED 

DYA 

= 

O.OOOOE+OO 

N-S/RAD 

CROSS-COUPLED 

DYB 

= 

O.OOOOE+OO 

N-S/RAD 

CROSS-COUPLED 

DAX 

= 

O.OOOOE+OO 

N-M-S/M 

CROSS-COUPLED 

DAY 

= 

O.OOOOE+OO 

N-M-S/M 

PRINCIPAL A 

DAA 

= 

O.OOOOE+OO 

N-M-S/RAD 

CROSS-COUPLED 

DAB 

= 

O.OOOOE+OO 

N-M-S/RAD 

CROSS-COUPLED 

DBX 

= 

O.OOOOE+OO 

N-M-S/M 

CROSS-COUPLED 

DBY 

- 

O.OOOOE+OO 

N-M-S/M 

CROSS-COUPLED 

DBA 

= 

O.OOOOE+OO 

N-M-S/RAD 

PRINCIPAL B 

DBB 

= 

O.OOOOE+OO 

N-M-S/RAD 
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GCYLT 


FILM THICKNESS DISTRIBUTION 

SAMPLE CASE 3: 3-LOBE GAS SEAL 


DIAMETER = 0-051 

M 

SPEED = 50000-00 RPM 

LENGTH = 0-025 

M 


CLEARANCE = 0-000013 

M 
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Figure 24. Clearance Distribution , Three-Lobe Seal 
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Figure 25. Pressure Distribution , Three-Lobe Seal 
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GCYLT 


PRESSURE 

DISTRIBUTION 

SAMPLE CASE 3: 3-LOBE GAS SEAL 



DIAMETER = 0.051 

M 

SPEED = 

50000.00 RPM 

LENGTH = 0-025 

M 



CLEARANCE = 0-000013 

M 
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Figure 26. Pressure Distribution, Viewing in Axial Direction 
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3.4 Sample Problem 4 - T-Shaped Sectored Seal 

This problem deals with an actual helium buffered seal analysis and design (for SSME) that was 
accomplished for NASA. A design that incorporates a self-adjusting clearance that can 
accommodate thermal and centrifugal distortions and shaft dynamic excursions avoids many of 
the problems associated with captured clearance configurations. The sectored ring seal provides 
the desired self-adjusting clearance features. The general configuration of the sectored seal is 
shown on Figure 27. The sectors consist of T-shaped sections mated to each other at each end 
with sealed joints. The sectors can move relative to each other circumferentially and that is how 
the seal accommodates variations in the sleeve dimensions due to thermal expansions and 
contractions and centrifugal growths. The T-shaped sector was chosen because it is a 
symmetrical shape, and the various fluid and friction forces can be designed to avoid upsetting 
moments on the individual sectors. An overlapping V joint prevents a direct clearance path 
between the hydrogen and oxygen ends of the seal. Each sector is supported by a hydrostatic 
fluid-film on its inner circumference and along the side walls which forms a friction-free 
secondary seal to permit free radial motion of the sectors in response to sleeve movements. The 
fluid films are predominantly hydrostatic to avoid any pitching tendencies introduced by the 
hydrodynamic effects. The hydrostatic bearings are fed by the buffer pressure on the outside 
diameter of the seal. Figure 28 shows the pressure distribution and force balance on the 
individual sectors. 

This sample problem describes one case conducted in the analysis of the circumferential 
hydrostatic seal on one of the sectors. The geometry and operating conditions are as follows: 

• Number of pads to be analyzed = 1 

• The position of the sector to satisfy a given load was to be determined 

• Load applied = 370 lb 

• Load angle from the x-axis = 270° 

• Initial guess on the eccentricity of the seal = 0.5 

• Initial guess on the eccentricity angle = 90° 

• Variable grids were used in both the axial and circumferential directions. The grid is 
made very fine around the source points. The starting angle of the sector is 30° from the 
x-axis and its angular extent is 120°. The axial length of the seal is 1.627 in. 

• Shaft diameter = 2.6798 in. 

• Reference clearance = 0.001 in. 

• Ratio of specific heat of the gas = 1.66 

• Gas constant = 1,790,000 in:/(s 2 -°R) 

• Absolute temperature = 528°R 

• Gas viscosity = 2.9 x 10 9 lb-s/in 2 

• Shaft speed = 0 rpm 

• Reference pressure = 14.7 psig 

• Boundary pressure surrounding the seal = 50 psig 

• Cross-coupled stiffness and damping were to be computed at an excitation frequency of 0 rpm 
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• There are 6 discrete, inherently compensated source points in the sector of diameter = 
0.020 in. The location of these orifices was determined from the design layout of the 
sector. The coefficient of discharge of each orifice is unity. 

• Buffer pressure = 200 psig 

• Flow is to be determined along four paths that make up the periphery of the seal. 

As indicated on the summary output (Table 4), the eccentricity of the sector to support the load 
is 0.6345 and the eccentricity angle is 90°. Plots of the clearance and pressure distribution are 
shown on Figures 29 and 30, respectively. Note on the plots the fine grid work surrounding the 
orifice locations. 



Figure 27. T-Shaped Sectored Ring Seal 
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F c = P s A + P a A a 

Radial Force Balance 
F r - F c ± F, ± F r = 0 

Axial Force Balance 

F ax - F p ' k s 8 a = 0 

871511-1 

Figure 28. Pressure Distribution and Force Balance, Sectored Seal 
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Table 4. Sectored Buffer Seal Sample Problem 


JOURNAL & LOAD POSITION 


ECCENTRICITY 


= 

ECCENTRICITY ANGLE 

= 

MINIMUM FILM 


= 

LOAD 


= 

LOAD ANGLE 


= 

POWER LOSS 


= 

LEAKAGE AT I = 1 


= 

LEAKAGE AT I = M 


= 

STIFFNESS COEFFICIENTS 


PRINCIPAL X 

KXX 

= 

CROSS-COUPLED 

KXY 

= 

CROSS-COUPLED 

KXA 

= 

CROSS-COUPLED 

KXB 

= 

CROSS-COUPLED 

KYX 

= 

PRINCIPAL Y 

KYY 

= 

CROSS-COUPLED 

KYA 

= 

CROSS-COUPLED 

KYB 

= 

CROSS-COUPLED 

KAX 

= 

CROSS-COUPLED 

KAY 

= 

PRINCIPAL A 

KAA 

= 

CROSS-COUPLED 

KAB 

= 

CROSS-COUPLED 

KBX 

= 

CROSS-COUPLED 

KBY 

= 

CROSS-COUPLED 

KBA 


PRINCIPAL B 

KBB 

= 

DAMPING COEFFICIENTS 


PRINCIPAL X 

DXX 

= 

CROSS-COUPLED 

DXY 

= 

CROSS-COUPLED 

DXA 

= 

CROSS-COUPLED 

DXB 

= 

CROSS-COUPLED 

DYX 

= 

PRINCIPAL Y 

DYY 

= 

CROSS-COUPLED 

DYA 

= 

CROSS-COUPLED 

DYB 

= 

CROSS-COUPLED 

DAX 

= 

CROSS-COUPLED 

DAY 

= 

PRINCIPAL A 

DAA 

= 

CROSS-COUPLED 

DAB 

= 

CROSS-COUPLED 

DBX 

= 

CROSS-COUPLED 

DBY 

= 

CROSS-COUPLED 

DBA 

= 

PRINCIPAL B 

DBB 

= 

RIGHTING MOMENT 

ABOUT X-X 

MX 


ABOUT Y-Y 

MY 

= 


- FLOW THRU SPECIFIED GRID 
FROM 1 1 TO 27 1 


-FLOW THRU SPECIFIED GRID 
FROM 1 34 TO 27 34 


-FLOW THRU SPECIFIED GRID 
FROM 1 1 TO 1 34 


-FLOW THRU SPECIFIED GRID 
FROM 27 1 TO 27 34 


0.63453 
90.00 DEG 
0.0003663 IN 
370.0 LB 
-90.00 DEG 

O.OOOOE+OO HP 

0.25673E-04 LB/S 

0.25673E-04 LB/S 


0.4032E+05 LB/IN 
-0.2428E-07 LB/IN 
-0.2531E-07 LB/RAD 
0.2009E-02 LB/RAD 
-141.6 LB/IN 
0.2429E+05 LB/IN 
-391.7 LB/RAD 
-74.14 LB/RAD 
0.8476E-06 IN-LB/IN 
-0.4463E-02 IN-LB/IN 
0.5055E+05 IN-LB/RAD 
-0.5231E-06 IN-LB/RAD 
0 . 1525E-02 IN-LB/IN 
0.1 127E-07 IN-LB/IN 
-0.4153E-08 IN-LB/RAD 
0. 1099E+05 IN-LB/RAD 


15.89 LB-S/IN 
0.3707E-10 LB-S/IN 
0.3617E-11 LB-S/RAD 
0.5322E-06 LB-S/RAD 
0.2025E-01 LB-S/IN 
128.5 LB-S/IN 
0.6261E-01 LB-S/RAD 
0.1058E-01 LB-S/RAD 
0.9535E-10 IN-LB-S/IN 
0.1738E-05 IN-LB-S/IN 
7.329 IN-LB-S/RAD 
-0 . 5647E-09 IN-LB-S/RAD 
-0.4686E-06 IN-LB-S/IN 
- 0 . 5672E- 1 2 IN-LB-S/IN 
-0.2274E-1 1 IN-LB-S/RAD 
1.288 IN-LB-S/RAD 


0.2995E-05 LB-IN 
0 . 2342E - 1 4 LB-IN 


LINE 

FLOW= -0.4651E-04 LB/S 


LINE 

FLOW= 0.4651E-04 LB/S 


LINE 

F(_OW= -0.2567E-04 LB/S 


LINE 

FLOW= 0.2567E-04 LB/S 
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Figure 29. Clearance Distribution, Sectored Seal 
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Figure 30. Pressure Distribution , Sectored Seal 
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3.5 Sample Problem 5 - Rayleigh-Step, Floating-Ring Seal 

This example represents another buffer fluid seal that was designed for use in the SSME. The 
principle of operation of a hydrodynamic, lift-pad, floating-ring seal is illustrated on Figure 3 1 . The 
seal consists of two rings that are mounted back-to-back. The buffer fluid enters between the rings 
and forces the rings up against the stationary housing. The buffer fluid leaks into the clearance 
annulus between the shaft and the seal and prevents ingress of exterior fluid on either side of the 
floating-ring assembly. The rings are held in equilibrium by a number of forces, as shown on Figure 
31. F c is a pressure force from the inlet buffer fluid that forces the rings up against the housings. 
This pressure force is partially balanced on the housing sides of the rings by undercutting and 
exposing the housing sides of the rings to buffer pressure. The balance force is identified as F B . F H 
represents a hydrodynamic force that is generated by rotation between the shaft and ring. The net 
hydrodynamic force is zero when the shaft and rings are in the concentric position. However, when 
the ring becomes eccentric with respect to the shaft, a hydrodynamic force is built up that opposes 
the eccentricity. There is also a normal force, F N , acting on the ring at the contact area between the 
ring and the housing. In addition to the equilibrium forces mentioned above, there is a friction force, 
F f , between the seal ring and housing. Seal tracking requires that the hydrodynamic force is strong 
enough to overcome the friction force. 

Figure 32 shows the hydrodynamic geometry that is incorporated into the bore of the seal rings. 
A portion of the length of the bore is segregated into sectors, and these sectors are separated from 
one another by axial grooves. A circumferential groove that goes completely around the bore 
is installed upstream of the final seal dam region. At the interior of the sectors, Rayleigh-step 
pockets are machined. The velocity direction of the shaft is such that it produces hydrodynamic 
pressures due to pumping of the fluid over the Rayleigh-step. The sealing occurs across the dam 
which is a narrow annulus of low clearance exposed to high pressure at its interior 
circumferential groove and to lower pressure at its outboard end. The shaded regions on Figure 
32 indicate depressions from grooves and Rayleigh-step pockets. 

In this example, one pad of the Rayleigh-step interface was examined from the high-pressure 
interior end to the low-pressure exterior end. Geometric and operating parameters are as follows: 

• International units are to be used 

• Because the performance of one pad was examined, the grid represented a 90° arc from 
the center of one axial groove to the next. The boundary conditions at the circumferential 
ends of the single pad must be periodic, i.e., all pads will act identically, which will occur 
when the shaft is in the concentric position. The periodicity option of the code was 
invoked. 

• Shaft diameter = 0.05 m 

• Total seal length = 0.0123 m 

• Seal clearance = 1.27 x 10~ 5 m 

• Step height = 2.54 x 10' 5 m 

• Gas viscosity = 2.19 x 10 s N-s/irf 

• Absolute temperature = 338. 6°K 

• Ratio of specific heat = 1.66 
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• Gas constant = 1154.8364 m 2 m 2 /(s 2 -°K) 

• Shaft speed = 70,000 rpm 

• Reference pressure = 101,352.93 Pa 

• The high pressure to be sealed is 1.37895 x 10 6 Pa, which would be at the bottom of the 
grid. The remaining boundaries are at 0 psig. 

In this example, advantage was taken of inputting constant pressure regions for the inlet grooves 
of sequential pads and the circumferential exterior grooves. 

Summary output for the single pad is shown on Table 5. Each individual pad load is 571.6 N. 
The leakage at the low-pressure end (I = M) exceeds that at the high-pressure end 1 = 1, because 
of the groove leakage. 

Following the output are the 3-D plots of the clearance and pressure distribution (Figures 33 and 
34, respectively). Note on Figure 33 that the high ambient pressures in the groove regions 
overshadow the contribution from the Rayleigh-step. 


Helium Buffer Fluid 



Figure 31. Rayleigh-Step , Floating Ring Seal 
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— Circumferential 
Groove 



Figure 32. 
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■Groove Depth = 0.229 mm 


Axial 

Groove 


^-Pocket- 7 

{Depth 

[ = 0.0254 mm) 


Section A-A 


50-mm Rayleigh-Step Pad 


Table 5. Single-Pad Rayleigh-Step Seal Problem 


-JOURNAL & LOAD POSI 
ECCENTRICITY 
ECCENTRICITY ANGLE 
MINIMUM FILM 
LOAD 

LOAD ANGLE 
POWER LOSS 
LEAKAGE AT I = 1 


ION 

= 0.00000 

= 0.00 DEG 

= 0.0000127 M 

= 571.6 N 

= -134.73 DEG 

= 26.27 W 

=-0.23044E-04 KG/S 


LEAKAGE AT I = M 

-RIGHTING MOMENT 
ABOUT X-X 
ABOUT Y-Y 


= 0.33187E-03 KG/S 

MX = -0.1609 N-M 

MY = 0.1552 N-M 
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Figure 33. Clearance Distribution, Rayleigh-Step Pad 
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Figure 34. Pressure Distribution , Rayleigh-Step Pad 
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3.6 Sample Problem 6 - Rayleigh-Step Seal with Eccentricity 

This problem will be similar to Problem 5 except the shaft is to be eccentric with respect to the 
seal ring. In this case, periodic boundary conditions cannot be used and to conserve grid space, 
one hydrodynamic pad will be modeled and the number of pads will be four. To model separate 
pads, however, requires that the boundary conditions be known on all extremities of the pad. The 
seal dam region is not a separate pad problem but is a single 360° pad. Thus, the problem 
resolves into two separate problems; one that treats the separate Rayleigh pads and one that treats 
the seal dam. For this particular example, only the Rayleigh-Step hydrodynamic region is 
considered. The following geometric and operating parameters have been applied: 

• International units were invoked 

• The shaft position relative to the seal ring was specified 

• Stiffness is calculated in four degrees of freedom at a 70,000 rpm excitation frequency 

• Number of pads is 4; each pad has an extent of 90° 

• Shaft diameter = 0.05 m 

• Shaft length = 0.0123 m 

• Reference clearance = 1.27 x 10 5 m 

• Gas viscosity = 2.19 x 10 5 N-s/m 2 

• Absolute temperature = 338.6°K 

• Ratio of specific heat = 1 .66 

• Gas constant = 1154.84 m 2 /(s 2 -°K) 

• Shaft eccentricity ratio = 0.5 

• Eccentricity angle = 270° 

• Shaft speed = 70,000 rpm 

• Reference pressure = 1.01353 x 10 5 Pa 

• Pad boundary pressures are 1.37 8 95 x 10 6 Pa. 

To include the groove pressures, two constant pressure regions were specified. 

As shown on Table 6, at the specified position, the load capacity of the seal is 69.73 N and the 
load angle is 69.47° from the x-axis. The minimum film thickness is 6.4 x 10 6 m. The clearance 
and pressure distributions are shown on Figures 35 and 36, respectively. 
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Table 6. Multipad Rayleigh-Step Seal Problem 

- JOURNAL & LOAD 

POSITION 


ECCENTRICITY 


= 0.50000 


ECCENTRICITY ANGLE 

= -90.00 

DEG 

MINIMUM FILM 


= 0.0000064 

M 

LOAD 


= 69.73 

N 

LOAD ANGLE 


= 69.47 

DEG 

POWER LOSS 


= 116.0 

W 

LEAKAGE AT I = 

1 

=-0.99207E-04 

KG/S 

LEAKAGE AT I = 

M 

= 0.99207E-04 

KG/S 

-STIFFNESS COEFFICIENTS 



PRINCIPAL X 

KXX 

= 0 . 1 298E+08 

N/M 

CROSS-COUPLED 

KXY 

= 0 . 5795E+07 

N/M 

CROSS-COUPLED 

KXA 

= 25.55 

N/RAD 

CROSS-COUPLED 

KXB 

= 2.247 

N/RAD 

CROSS-COUPLED 

KYX 

= 0.2313E+06 

N/M 

PRINCIPAL Y 

KYY 

= 0. 1616E+08 

N/M 

CROSS-COUPLED 

KYA 

= 59.29 

N/RAD 

CROSS-COUPLED 

KYB 

= 30.28 

N/RAD 

CROSS-COUPLED 

KAX 

= 0.2214E-07 

N-M/M 

CROSS-COUPLED 

KAY 

= -0.9859E-08 

N-M/M 

PRINCIPAL A 

KAA 

= 36.49 

N-M/RAD 

CROSS-COUPLED 

KAB 

= 0.2037 

N-M/RAD 

CROSS-COUPLED 

KBX 

= 0. 1360E-07 

N-M/M 

CROSS-COUPLED 

KBY 

= -0.5505E-07 

N-M/M 

CROSS-COUPLED 

KBA 

= -12.66 

N-M/RAD 

PRINCIPAL B 

KBB 

= 19.24 

N-M/RAD 

-DAMPING COEFFICIENTS 



PRINCIPAL X 

DXX 

= 753.6 

N-S/M 

CROSS-COUPLED 

DXY 

= -303.8 

N-S/M 

CROSS-COUPLED 

DXA 

= -0. 1023E-02 

N-S/M 

CROSS-COUPLED 

DXB 

= -0.5204E-04 

N-S/RAD 

CROSS-COUPLED 

DYX 

= 163.4 

N-S/M 

PRINCIPAL Y 

DYY 

= 815.2 

N-S/M 

CROSS-COUPLED 

DYA 

= -0.2082E-02 

N-S/RAD 

CROSS-COUPLED 

DYB 

= -0. 1277E-02 

N-S/RAD 

CROSS-COUPLED 

DAX 

= -0.2680E- 14 

N-M-S/M 

CROSS-COUPLED 

DAY 

= 0.3079E- 13 

N-M-S/M 

PRINCIPAL A 

DAA 

= 0.3602E-02 

N-M-S/RAD 

CROSS-COUPLED 

DAB 

= -0. 1509E-03 

N-M-S/RAD 

CROSS-COUPLED 

DBX 

= -0.2348E-13 

N-M-S/M 

CROSS-COUPLED 

DBY 

= 0.8120E- 13 

N-M-S/M 

CROSS-COUPLED 

DBA 

= -0.51 16E-04 

N-M-S/RAD 

PRINCIPAL B 

DBB 

= 0.2658E-02 

N-M-S/RAD 

-RIGHTING MOMENT 




ABOUT X-X 

MX 

= -0. 1978E- 15 

N-M 

ABOUT Y-Y 

MY 

= 0.6593E- 16 

N-M 
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GCYLT 

RRYLEIGH - STEP SERE PROBLEM 
DIRMETER = 0-050 M 

LENGTH = 0-012 M 

CLERRRNCE = 0-000013 M 


FILM THICKNESS DISTRIBUTION 


SPEED = 70000-00 RPM 



Figure 35. Clearance Distribution , Rayleigh-Step Seal with Eccentricity 
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95TR34-V2 


Figure 36. Pressure Distribution, Rayleigh-Step Seal with Eccentricity 
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3.7 Critical Mass - Sample Problem 

The Rayleigh-Step Seal problem (see Sections 3.5 and 3.6) was further analyzed to determine 
critical mass and frequency as a function of speed and operating clearance. 

Figure 37 shows critical mass as a function of speed and operating clearance. If the mass 
attributable to the seal or bearing exceeds the critical mass, then an instability can occur. The 
usefulness of the critical mass parameter is that it provides a comparative measure of the stability 
characteristics of different configurations. Clearly from Figure 37, the low clearance seal would 
have superior stability characteristics. The critical frequencies are shown on Figure 38. For both 
clearances, a constant ratio exists between the critical orbital frequency and the operating speed. 
The ratios are slightly less than 0.5. 

Critical mass problems should be confined to small axial pressure gradient cases, because the 
angular modes are not considered. 



Figure 37. Rayleigh-Step Seal, Critical Mass versus Speed 
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Figure 38. Rayleigh-Step Seal, Critical Frequency versus Speed 
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4.0 COMPARISONS OF GCYLT WITH GCYL 


It is informative to determine the effects of turbulence on seal performance. Comparisons were 
made between the laminar code GCYL and the revised turbulent code GCYLT to establish the 
significance of turbulence. The sample problems previously described were run for each of the 
codes and the results are shown in Tables 7 through 12. 

Table 7 shows the variation for the single Rayleigh-Step problem described by Sample Problem 1 . 
The Couette Reynolds number for this case is 


, PVh 

p 0 RaC 

e V 

GT\i 

c a~ 

( 1.9685] 

{ 70,000* 2tc ) 

i 2 J 

l 60 J 


(250,000) (530) (3a 10" 9 ) 


(. 001 ) 


= 3896 


Thus, the problem is clearly in the turbulent regime (R e >1000) and accounts for the variation in 
performance. In particular, note the higher load capacity of the turbulent model. The flow is 
also greater because of the higher pressure generated. 


Table 7. Rayleigh-Step Seal Comparison 



GCYLT 

GCYL 

£ 

Y> deg 

h min. in. 

W, lb 
a L deg 
HP, hp 
Q„ lb/s 
Lnax’ lb/in? 

.40 

-90 

0.0006015 

131.5 

58.49 

2.308 

-0.10191 x 10' 2 
54.325 

0.4 

-90 

0.0006015 

27.54 

61.44 

0.4555 

-0.44783 x 10“ 3 
12.209 


£ = 

eccentricity ratio 

y 

eccentricity angle 

hmin 

minimum film thickness 

w 

load 

<*L = 

load angle 

HP = 

power loss 

Q. = 

leakage at I =1 

p m; , = 

maximum film pressure 
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Table 8 shows the comparison for the three-pad sectored seal. In this case, the Couette Reynolds 
number of 158 is definitely laminar. However, the externally pressurized orifices introduce high 
pressure gradients which cause Poiseuille turbulence to occur and cause a slight variation in results. 

The three-lobe seal comparisons are shown on Table 9. The Couette Reynolds number is 841, 
which is in the transition zone and which will cause a slight variation with the laminar results. 

The hydrostatic sector comparison (Sample Problem 4) is shown on Table 10. Since this is a 
zero speed case, Couette turbulence could not enter. Poiseuille turbulence, however, is present 
and was confirmed by a temporary flag inserted into the code. The turbulence occurs at the 
orifice locations, where high pressure gradients occur and cause the variation in performance 
outlined in the tabulation. 

The last two examples (Tables 11, 12) were for Rayleigh-step seals operating under similar 
conditions. The Couette Reynolds number for the step region is 384 and is, thus, laminar. The 
variation between the laminar and turbulent cases is due to Poiseuille turbulence occurring in the 
step region where high pressure gradients occur. 


Table 8. Three-Pad Sectored Seal Comparison 



e = eccentricity ratio 

y ~ eccentricity angle 

h min = minimum film thickness 

W = load 

a L = load angle 

HP = power loss 

Q, = leakage at I- 1 

P max = maximum film pressure 
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Table 9. Three-Lobe Seal Comparison 



GCYLT 

GCYL 

8 

.22088 

.22103 

y, deg 

129.42 

129.38 

h m,n> in- 

0.0000037 

0.0000037 

W, lb 

200.2 

200.2 

a h , deg 

-90 

-90 

HP, hp 

186.2 

186.1 

Q„ lb/s 

.12898 x lO' 4 

.12875 x 10 4 

K xx , lb/in. 

.152 x 10 9 

.146 x 10 9 

K xy , lb/in. 

-.2103 x 10 8 

-.2379 x 10 8 

K yx , lb/in. 

-.3595 x 10 8 

-.3890 x 10 8 

Kyy, lb/in. 

.8662 x 10 8 

.1002 x 10 9 

D xx , lb-s/in. 

.1165 x 10 5 

9939 

D xy , lb-s/in. 

-.1510 x 10 5 

-.1089 x 10 5 

D yx , lb-s/in. 

5204 

5140 

Dyy, lb-s/in. 

.1653 x 10 5 

.1393 x 10 5 


Kjj = Stiffness in i direction due to j displacement 

Djj = Damping in i direction due to j velocity 

e = eccentricity ratio 

Y = eccentricity angle 

h mjn = minimum film thickness 

W - load 

a L = load angle 

HP = power loss 

Q[ = leakage at 1=1 

P maj[ = maximum film pressure 
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Table 10. Hydrostatic Sector Comparison 



GCYLT 

GCYL 

e 

0.6345 

0.55932 

y, deg 

90 

90 

hmm> in. 

0.000366 

0.0004414 

W, lb 

370 

369.6 

a L , deg 

-90 

-90 

HP, hp 

0 

0 

Q„ lb-s 

-0.25673 x 10' 4 

-0.39586 x 10' 4 

Qm. lb ' s 

0.25673 x 10' 4 

0.39586 x 10' 4 

K, x , lb/in. 

.40320 

43020 

Ky X , lb/in. 

-141.6 

-126.9 

Kyy, lb/in. 

24290 

48580 

Ky a , lb/rad 

-391.7 

-298.6 

Ky b , Ib/rad 

-74.14 

-59.98 

K AA , in. -lb/rad 

50550 

43,890 

K bb , in. -lb/rad 

10990 

10,000 

D xx , lb-s/in. 

15.89 

10.98 

D yx , lb-s/in. 

0.0202 

0.0143 

Dyy, lb-s/in. 

128.5 

81.68 

D ya , lb-s/rad 

0.063 

0.0367 

D yb , lb-s/rad 

0.0106 

0.0067 

D aa , in. -lb-s/rad 

7.329 

4.599 

D bb , in.-lb-s/rad 

1.288 

0.9138 

P„ lb/in? 

178 

170* 

P 2 , lb/in? 

193 

188* 

P„ lb/in? 

197 

194* 

P 4 , lb/in 2 

197 

194* 

P 5 , lb/in 2 

193 

188* 

P 6 , lb/in 2 

178 

170* 


^Pressures downstream of orifices 

e = eccentricity ratio 

a L = load angle 

y = eccentricity angle 

HP = power loss 

h mi[| = minimum film thickness 

Q, = leakage at 1=1 

W = load 

Q v , = leakage at I=M 

P max = maximum film pressure 
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Table 11. Rayleigh-Step Seal , Single Pad Comparison 



GCYLT 

GCYL 

e 

y, deg 

K>v m 
W, N 
a L , deg. 
HP, hp 
Q„ kg/s 
Qm> kg/s 

0 

0 

1.27 x 10 5 
571.6 
-134.73 

26.27 

-0.23 x lO' 4 
.332 x 10 3 

0 

0 

1.27 x lO 5 
571 

-134.78 

21.79 

-.217 x 10‘ 4 
.403 x 10' 3 


e = 

eccentricity ratio 

1 

eccentricity angle 

Kin = 

minimum film thickness 

w = 

load 

«L = 

load angle 

HP = 

power loss 

Q, = 

leakage at 1=1 

Pmin = 

maximum film pressure 
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Table 12. Four-Pad Rayleigh-Step Seal Comparison 



GCYLT 

GCYL 

e 

0.5 

0.5 

7, deg 

-90 

-90 

h min - in. 

6.4 x 10' 6 

6.4 x 10‘ 6 

W, N 

69.73 

52.96 

a L , deg 

69.47 

71.63 

HP, hp 

116 

95.08 

Qi, kg/s 

-.99207 x 10 4 

-.88483 x lO' 4 

Qm- kg/s 

.99207 x 10' 4 

.88483 x 10" 4 

K xx , N/in. 

12,980,000 

10,160,000 

K xy , N/in. 

5,795,000 

3,371,000 

K xa , N/rad 

25.55 

14.33 

K xb , N/rad 

2.47 

-.3405 

K yx , N/in. 

23.130 

-61,600 

Ky y , N/in. 

16,160 

11,690 

K ya , N/rad 

59.29 

37.51 

K yb , N/rad 

30.28 

20.59 

K aa , N-m/rad 

36.49 

27.53 

K AB , N-m/rad 

0.2307 

1.657 

K ba , N-m/rad 

-12.66 

-10.32 

K bb , N-m/rad 

19.24 

14.46 

D xx , N-s/in. 

753.6 

786.4 

D xy , N-s/in. 

-303.8 

-216.6 

D yx , N-s/in. 

163.4 

192.7 

Dyy, N-s/in 

815.2 

901.7 


e = eccentricity ratio 

a L - load angle 

y = eccentricity angle 

HP - power loss 

h min = minimum film thickness 

Q, = leakage at 1=1 

W = load 

P max = maximum film pressure 
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5.0 VERIFICATION 


Several mechanisms were used to conduct verification of the code. Results of the code were 
compared against information in the public domain literature, and comparisons were made against 
the results of other codes and against manual computations. 

Although simple in concept, extensive changes were implemented to the original laminar code 
to accommodate turbulence. A first check on the turbulent code was to run a laminar case to see 
if it compared precisely against the original code. 

Most of the comparison cases involved a 360° plain seal with the parameters shown on Table 13. 


Table 13. 360° Plain Seal Parameters 


Diameter 

Length 

Clearance 

Viscosity 

Absolute Temperature 
Specific Heat Ratio 
Gas Constant 
Eccentricity Ratio 
Eccentricity Angle 
Ambient Pressure 
Pressure Difference 
Speed 
W 
Y 

HP 

K* 

D:; 


1 in. 

2 in. 

0.0005 in. 

3 x 10' 9 lb- s/in 2 
530°R 

1.4 

247,000 in 2 /S 2 /°R 

0.5 

270° 

Pa = Variable (psia) 

Pd = Variable (psi) 

N = Variable (rpm) 
Load Capacity (lb) 
Attitude Angle (deg) 
Power Loss (hp) 
Stiffness in i direction 
due to j displacement 
Damping in i direction 
due to j velocity 


Laminar comparisons were made against the original codes and against another code, SPIRALG, 
that was run with the smooth surface option. Results are shown on Table 14. The turbulent code 
GCYLT compared exactly against GCYL and comparisons with SPIRALG were excellent. 

The original code GCYL was also compared against information in the literature. Figure 39 shows 
the pressure distribution in an infinitely long Rayleigh-Step slider at various values of the 
compressibility parameter A., For this case, closed-form solutions are available. Turbulence would 
show some variation at high values of A. 
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Table 14. Laminar Check of GCYLT 



GCYL 

GCYLT 

SPIRALG 

^min’ 

mils 

0.000252 

0.000252 

.00025 

w, 

lb 

1.695 

1.695 

1.672 

Y 

deg 

7.76 

7.76 

7.47 

HP 

hp 

.214 x 10 4 

.214 x 10‘ 4 

.214 x 10‘ 4 

K xx 

lb/in. 

921 

921 

869 

K xy 

lb/in. 

8572 

8572 

8437 

Ka 

lb/rad 

2.163 

2.163 

— 

KxJ3 

lb/rad 

2.431 

2.431 

— 

Kyx 

lb/in. 

-6715 

-6715 

-6630 

Kyy 

lb/in. 

2301 

2301 

2166 

Kya 

lb/rad 

1.072 

1.072 

— 

k, b 

lb/rad 

0.5987 

0.5987 

— 

K ax 

in.-lb/in. 

— 

— 

-- 

If 

AV ay 

in.-lb/in. 

— 

— 

— 

K aa 

in. -lb/rad 

329.4 

329.4 

293 

K aB 

in. -lb/rad 

819.3 

819.3 

804 

Kb x 

in.-lb/in. 

— 

— 

— 

Kfiy 

in.-lb/in. 

— 

— 

- 


in.-lb/rad 

-1399 

-1399 

1369 

k bb 

in.-lb/rad 

97.87 

97.87 

86.33 

D xx 

lb-s/in. 

125.2 

125.2 

124 

Dxy 

lb-s/in. 

-30.16 

-30.16 

-29 

D X a 

lb-s/rad 

— 

-- 

-- 

^xB 

lb-s/rad 

— 

— 

— 

D yx 

lb-s/in. 

34.73 

34.73 

33.87 


lb-s/in. 

205.1 

205.1 

202.9 

Dya 

lb-s/rad 

— 

— 

— 

DyB 

lb-s/rad 

— 

— 

— 

D ax 

in 2 -lb/s 

— 

— 

-- 

D ay 

in 2 -lb/s 

— 

— 

— 


in.-lb-s/rad 

30.76 

30.76 

30.31 

^ai! 

in.-lb-s/rad 

-2.787 

-2.787 

-2.58 

D B x 

in 2 -lb-s 

— 

— 

— 

D By 

in 2 -lb-s 

— 

— 

— 


in.-lb-s/rad 

2.698 

2.698 

2.48 

^BB 

in.-lb-s/rad 

16.99 

16.99 

16.75 


P a = 14.7 psia P d = 0.0 N = 1000 rpm Excitation Frequency = 0 
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Ll 


h 2 

• GCYL I f 

— Exact Solution* 

A = 6y UL/p a h \ 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


x/L 


86259-1 


* "Theory of Hydrodynamic Lubrication", O. Pinkus, R. Stemlicht, McGraw-Hill, N.Y. 1961 


Figure 39. Rayleigh- Step, Program Verification 
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Further comparisons were made for a plain cylindrical seal with an L/D ratio of 1 with 
information from Reference 8. Computations were made at two different eccentricity ratios, e 
= 0.6 and 0.8. Nondimensional load capacity and attitude angles are shown in Figures 40 and 
41, respectively. Excellent correlation is demonstrated. Again, these comparisons were made 
with the earlier laminar version of the code. 

To validate Couette and Poiseuille turbulence, GCYLT was run against another code available 
at MTI called GBEAR. The GBEAR code accurately predicts Couette turbulence for gases, but 
not Poiseuille turbulence. 

Table 15 shows comparative results for Couette turbulence with a Reynolds number of 48,995. 
The comparative results are very good for all parameters. 

A Poiseuille turbulence check was made by running at a low speed, with a high pressure gradient 
and a high ambient pressure. By using a high value of ambient pressure, the compressibility 
parameter, A, is reduced to a low value. Then the liquid and gas cases nearly correspond. 
Comparative results are shown on Table 16. Again, the comparative results are very good. 


Plain Cylindrical Seal, L/D = 1 



Figure 40. Dimensionless Load Capacity versus A 
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A = 6|Uo(R/C)**2/P 

Figure 41. Attitude Angle versus A 


10 
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Table 15. Couette Turbulence Comparison 



GCYLT 

GBEAR 

L 

A1 min 

mils 

0.252 

0.252 

w 

lb 

1270 

1269 

7 

deg 

4.92 

4.63 

HP 

hp 

1.2 

1.2 

Q 

Ib/s 

- 

— 

K» 

lb/in. 

43,740 

40,917 

K xy 

lb/in. 

6,458,000 

6,443,64 

K 

lb/in. 

-5,062,000 

5,058,81 

Kyy 

lb/in. 

919,000 

813,689 

D,. 

lb-s/in. 

1920 

1936 

D xy 

lb-s/in. 

-270.7 

-166 

Dyx 

lb-s/in. 

334.3 

-147 

Dyy 

lb-s/in. 

2649 

2650 

M c 

lb 

16 

12 

C0 c 

rpm 

23,987 

24,064 


N = 50,000 rpm P a = 14,700 psig 

Rec = 48,995 Excitation Frequency = 0 


Table 16. Poiseuille Turbulence Comparison 



GCYLT 

GBEAR 

L 

min 

mils 

0.25 

0.25 

w 

lb 

9.513 

9.449 

7 

deg 

0.03 

0 

HP 

hp 

9.48 x 10' 5 

9.43 x 10 5 

Q. 

lb/s 

.0270 

.0267 

K xx 

lb/in. 

10.05 

— 

K xy 

lb/in. 

46,170 

47,271 

Kyx 

lb/in. 

-38,860 

-37,188 

Kyy 

lb/in. 

23.16 

0.78 

D xx 

lb-s/in. 

728.7 

739 

D xy 

lb-s/in. 

-0.60 

- 

D yx 

lb-s/in. 

0.82 

0.88 

D yy 

lb-s/in. 

862.7 

875 

M c 

lb 

-.35 

-2.8 

C0 c 

rpm 

511 

497 


N = 1000 rpm P a = 14,700 psig 

P d = 500 psig Excitation Frequency = 0 
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An internal check of the code was made by analyzing a recessed hydrostatic bearing. With the 
flow path option, the net flow around the periphery of a hydrostatic pad can be determined and 
compared against the inflow to the recess. For flow continuity, the sum of the peripheral flows 
should equal the inlet flow. The following geometry and operating parameters were considered. 

• A single pad with grid dimension of 15 x 37 (M x N) 

• The pad diameter is 2 inches 

• The pad length is 2 inches 

• The pad clearance is 0.001 in. 

• The pad angle is 180° and the starting angle is at 180° 

• There is one recess located in the pad, and the grid comer points are as follows: 

Left bottom corner, M = 3, N = 22 

Right top comer, M = 13, N = 27 

• The specific heat of the gas is 1 .4 

• The gas constant is 250,000 in?/(s 2 -°R) 

• The absolute temperature is 530°R 

• The absolute viscosity is 3 x 10" 9 Ib-s/in? 

• The inlet orifice diameter to the recess is 0.020 in. and the coefficient of discharge is 1.0. 
The orifice is located in the grid at M = 8, N = 24. 

• The supply pressure to the orifice is 150 psig. The pressure surrounding the pad is at 0 
psig. The reference ambient pressure is 14.7 psia. 

• Several eccentricities and speeds were examined and are defined in the subsequent 
discussions. 

The output from the code supplies the total flow from the peripheral flow path and the pressure 
in the recess. A manual computation can then be made for calculating the inlet flow through the 
orifice using the following equation: 


4 


386.4 A C n G t P 

o D 1 s 


r 

r "y 

2 


f y 
Pj_ 

Y-l 



Pj_ 

Y 

1- 

Y 



S'. 



S'. 




(5-1) 


where 


G. = 

2y 


G 0(y -1) 

f 0 = inlet flow, Ib/s 

A 0 = orifice area, in? 

C D = discharge coefficient 

p s = supply pressure, psia 

p r = recess pressure, psia 

Y = ratio of specific heats 

G c - gas constant, in:/(s 2 -°R) 

0 = absolute temperature. 
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Table 17 provides the results of several cases: 


Table 17. Recessed Pad Flow Comparisons 



N 

Q P 

P r 

Qo 

A 

e 

rpm 

lb/s 

psig 

Ib/s 

% 

0.0 

0.0 

0.001188 

46.6* 

0.001189 

0.08 


0.0 

0.001162 

88.6 

0.001163 

0.09 

0.0 

70,000 

0.001188 

37.7* 

0.001189 

0.08 


* = choked flow 

N = shaft speed 
P r = recess pressure 
A = percent variation 


e = eccentricity ratio 
Q p = peripheral flow 
Q 0 = orifice flow 


Note that the peripheral and orifice flows differ by less than 0.1%. 

When using the source points or spot recess options of the code, it is important to surround the 
source point with a fine grid to obtain an accurate result and a computation in which pressures 
will converge. Studies were made of varying grid sizes for a source problem. The variable grid 
option was applied and varied. A single pad with a central row of orifices were analyzed (see 
Sample Problem Number 4). The following information is pertinent: 

• Number of pads = 1 

• Pad angle =120° 

• Start angle = 30° 

• Number of grid points in circumferential direction = 34 

• Number of grid points in axial direction = 27 

• Diameter = 2.6798 in. 

• Length = 1.627 in. 

• Specific heat ratio = 1.66 

• Gas constant = 1 ,790,000 im/(s 2 -°R) 

• Absolute temperature = 528 °R 

• Viscosity = 2.9 x 10" 9 lb-s/in: 

• Shaft speed = 0 rpm 

• Reference pressure = 14.7 psia 

• Boundary pressure = 50 psig 

• Supply pressure to inherently compensated orifices = 200 psig 

• Preload = 50% located at the center of the pad 
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• Stiffness is to be determined 

• Six source points are located along a circumferential line in the axial center of the pad 
at circumferential grid locations 5, 10, 15, 20, 25, 30. 

The hole diameter is 0.020, and the coefficient of discharge is 1 .0. 

Table 18 indicates the effect of grid width around the source point in both the axial and 
circumferential directions. As the grid width is changed, the source pressures remain relatively 
unaffected until the grid width is 6 to 8X the orifice hole size. A similar conclusion can be drawn 
for the other performance parameters of load, flow, stiffness and damping. The recommended grid 
width from the source point to a neighboring grid line is twice the orifice diameter. 


Table 18. Comparative Studies - Discrete Orifices versus Grid Size A 
(Orifice Size = 0.015 in.) 


Comparison of Source Pressures 


A 

in 


P, 

psig 


1 2 
psig 


A 3 

psig 


1 4 

psig 


1 5 

psig 


1 6 

psig 


0.015 

0.030 

0.060 

0.120 


166 

170 

166 

160 


187 

185 

184 

182 


192 

191 

190 

188 


192 

191 

190 

189 


187 

185 

184 

182 


174 

170 

166 

160 


Comparison of Performance 


A 

in. 

W 

lb 

Q, 

lb/s 
x 10 4 

Qm 
lb/s 
x 10 4 

lb/in. 
x 10' 6 

Kyy 

Ib/in. 
x 10" 6 

D xx 

(lb-s)/in. 

D yy 

(lb-s)/in. 

0.015 

356.9 

0.48753 

0.48753 


0.0858 

8.506 

60.67 

0.030 

361.7 

0.50539 

0.50539 

1 

0.0842 

8.494 

59.79 

0.060 

368.2 

0.5305 

0.5305 

0.0475 

0.0817 

8.360 

58.70 

0.120 

377.3 

0.5651 

0.5651 

0.0522 

0.0967 

8.072 

55.97 


A = grid width in both circumferential and axial directions 

W = load capacity 

Q, = flow out of grid line M=1 

Q V1 = flow out of grid line M=M 

K xx and K = stiffness in x and y directions, respectively 

D xx and D = damping in x and y directions, respectively 
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6.0 THEORETICAL DESCRIPTION AND NUMERICAL 
METHODS FOR CODE GFACE 

6.1 Governing Equations 

Reynolds lubrication equation for compressible flow in polar coordinates is as follows: 


3 ( rph 3 dp} 1 dpf p h 3 3/?^ 
r 30 [ |i 30 


3 r 


V 


(i dr 


=6t/ M +12n ^ 

30 3 1 


(6-1) 


The coordinate system used is shown on Figure 42. The equation is made dimensionless with the 
following definitions. (Upper case variables are dimensionless). 


r 

R = —- H- 

h 

; T=-- P = -?-9 

P 

C 0 

K Pc 

A _ 6(4 ay; 2 

t ■- 

- 12^; 

p.cl ’ 

o 

PcCl 

Substituting 

the 

dimensionless variables 


Reynolds equation: 


(6-2) 


Reynolds equation produces a dimensionless 


3 

3 R 


f 

RPH 3 
V 



J__3_ 
R 30 



. d{PH) d(PH ) 

A 1 

30 3 T 


(6-3) 


For steady-state solutions, the time dependent term on the right-hand side is eliminated. 

The method of solutions is equivalent to that described for GCYLT. The grid is portioned into 
cells and a flow balance is established throughout the cells to determine the pressure distribution 
utilizing Newton-Raphson in combination with the column-matrix approach for each iteration. 
Precise details are included in Reference 3. 
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Figure 42 . Coordinate System for GFACE 
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6.2 Film Thickness Distribution and Misalignment 

Misalignment occurs through axes in the plane of the bearing. If a thrust collar misaligns an angle 
a about the x axis and P about the y axis (see Figure 42), the film thickness at any radius r, and 
angle 0 is : 


H- C Q + e - ar sin§ + $r cos§ 

where C Q - reference or initial clearance (6-4) 

e =axial displacement 


In non-dimensional form 


H= 1 + - — -^tfsine + ^^cose 

c. c. c. 


(6-5) 


In addition to a flat misaligned surface, other geometries can be treated. These are depicted on 
Figure 43. Included are Rayleigh steps and tapers in both the circumferential and radial 
directions. The boundaries of circumferential steps and tapers are defined by the ij grid point 
values of the lower left and upper right hand corners of the depressed region. Hydrostatic 
geometries are also included. Inherently compensated orifices can be located in the grid region as 
well as orifice fed holes and orifice fed recesses. The recess location is defined similarly to steps 
and tapers. 
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Section D-D 

• Circumferential Taper 




Section C-C 

• Rayleigh Step 

• Section E - E 

• Hydrostatic Recess 


87907-1 


Figure 43. Non-Uniform Clearance Geometries 
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7.0 SAMPLE PROBLEMS OF CODE GFACE 


7.1 Multipad Rayleigh Step Seal 

An eight-pad Rayleigh-step seal without pressure gradient, but with misalignment was analyzed. 
The geometry and operating conditions were as follows: 


Number of pads, NPAD 

Outer diameter 

Inner diameter 

Angular extent of pad 

Starting angle of first pad 

The displacement was given 

Compute Stiffness at synchronous frequency 

Apply variable grid 

Clearance 

Misalignment angle about the X- axis 

Number of steps in grid 

Step depth 

Location of step 

Lower left corner 

Upper right comer 

Specific heat ratio 

Gas constant 

Absolute temperature 

Viscosity 

Speed 

Convergence tolerance 
Ambient (reference) pressure 


8 

4.0 in. 

2.5 in. 

35 degrees 
5 degrees 

5000 

0.0002 in. 

-0.001 degrees 
1 

0.0004 in. 

1=5, J=1 
1=8, J=8 
1.4 

246,900 in 2 /s 2 /°R 
1460° R 

5.35 x 10- 9 lb- s/in 2 
5000 rpm 
0.01 
14.7 psi 


Summary output is indicated on Table 19. 

Clearance and pressure distribution plots are shown on Figures 44, 45 and 46. Since 
misalignment is slight, the variation from pad to pad is not very significant. 


7.2 Tapered Land Seal 

A single sector of a 16 pad tapered land seal was analyzed. Periodic boundary conditions were 
applied. Figure 47 shows the grid network layout of the pad. High pressure at the top flows 
through a radial groove into the tapered region and also into a circumferential groove that 
circumvents the complete seal. The slope of the taper is 0.05209. 
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• A load of 13 lb was applied and the axial position to support this load was determined. 
The initial guess for clearance was 0.0002 in. 

• Stiffness and damping coefficients were determined at zero speed excitation. 

• Periodic boundary conditions were applied to a single pad representation. The start of the 
pad angle was 67.5 degrees from the x-axis and the pad extent was 22.5 degrees. 

• The outer diameter was 4.5 in. and the inner diameter was 3.793 in. 

• As shown on Figure 47, the slope of the taper is 0.0509 and its location is indicated on 
the figure. The radial inlet region and circumferential groove were modeled as 
constant pressure areas. 

• The gas viscosity was 1.75 x 10" 9 lb-s/in. 2 and the gas constant was 423,184 in. 2 /s 2 /°F. 

• The operating speed was 50,000 rpm and the pressure differential was 50 psig from 
OD to ID. 

The program output is indicated on Table 20. The computed load is 13.08 lb as compared to the 
13 lb specified, which is well within the tolerance. The computed clearance is 0.000296 in. Plots 
of the clearance and pressure distribution are indicated on Figures 48 and 49. 


Table 19. GFACE Multipad Step Seal 


-JOURNAL POSITION & FORCES 


MINIMUM CLEARANCE 


= 

0.000165 

IN 

LOAD CAPACITY 


= 

39.43 

LB 

MOMENT ABOUT X-X 


= 

6.680 

IN-LB 

MOMENT ABOUT Y-Y 


= 

-1.503 

IN-LB 

POWER LOSS 


= 

0. 1527E-01 

HP 

LEAKAGE AT INNER 


= 

-0.6470E-06 

LB/S 

LEAKAGE AT OUTER 


= 

0.1003E-05 

LB/S 

STIFFNESS COEFFICIENTS 




PRINCIPAL Z 

KZZ 

= 

0.2805E+06 

LB/IN 

CROSS-COUPLED 

KZA 

= 

0.5135E+05 

LB/RAD 

CROSS-COUPLED 

KZB 

= 

-5180. 

LB/RAD 

CROSS-COUPLED 

KAZ 

= 

0.5199E+05 

IN-LB/IN 

PRINCIPAL A 

KAA 

= 

0.3825E+06 

IN-LB/RAD 

CROSS-COUPLED 

KAB 

= 

0.4896E+05 

IN-LB/RAD 

CROSS-COUPLED 

KBZ 

= 

-0. 1025E+05 

IN-LB/IN 

CROSS-COUPLED 

KBA 

= 

-0.5135E+05 

IN-LB/RAD 

PRINCIPAL B 

KBB 

= 

0.3771E+06 

IN-LB/RAD 

DAMPING COEFFICIENTS 




PRINCIPAL Z 

DZZ 

= 

17.97 

LB-S/IN 

CROSS-COUPLED 

DZA 

= 

-18.28 

LB-S/RAD 

CROSS-COUPLED 

DZB 

= 

-6.224 

LB-S/RAD 

CROSS-COUPLED 

DAZ 

= 

-19.62 

IN-LB-S/IN 

PRINCIPAL A 

DAA 

= 

19.59 

IN-LB-S/RAD 

CROSS- COUPLED 

DAB 

= 

-89.65 

IN-LB-S/RAD 

CROSS-COUPLED 

DBZ 

= 

0. 1520E-02 

IN-LB-S/IN 

CROSS-COUPLED 

DBA 

= 

90.14 

IN-LB-S/RAD 

PRINCIPAL 8 

DBB 

= 

24.25 

IN-LB-S/RAD 
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Table 20. GFACE Single Tapered Land Seal 

-JOURNAL POSITION & FORCES 


MINIMUM CLEARANCE 


= 

0.000296 

IN 

LOAD CAPACITY 


= 

13.08 

LB 

MOMENT ABOUT X-X 


= 

26.82 

IN-LB 

MOMENT ABOUT Y-Y 


= 

-5.263 

IN-LB 

POWER LOSS 


= 

0.2926E-01 

HP 

LEAKAGE AT INNER 



-0.2327E-04 

LB/S 

LEAKAGE AT OUTER 


= 

0.9632E-05 

LB/S 

STIFFNESS COEFFICIENTS 




PRINCIPAL Z 

KZZ 

- 

4719. 

LB/IN 

CROSS-COUPLED 

KZA 

= 

9396. 

LB/RAD 

CROSS-COUPLED 

KZB 

= 

-2108. 

LB/RAD 

CROSS-COUPLED 

KAZ 

= 

9939. 

IN-LB/IN 

PRINCIPAL A 

KAA 

= 

0. 1990E+05 

IN-LB/RAD 

CROSS-COUPLED 

KAB 

= 

-4320. 

IN- LB/RAD 

CROSS-COUPLED 

KBZ 

= 

-1544. 

IN-LB/IN 

CROSS-COUPLED 

KBA 

= 

-2980. 

IN-LB/RAD 

PRINCIPAL B 

KBB 

= 

1237. 

IN-LB/RAD 

DAMPING COEFFICIENTS 




PRINCIPAL Z 

DZZ 

= 

1.608 

L8-S/IN 

CROSS-COUPLED 

DZA 

= 

3.256 

LB-S/RAD 

CROSS-COUPLED 

DZB 

= 

-0.6525 

LB-S/RAD 

CROSS-COUPLED 

DAZ 

= 

3.255 

IN -LB-S/IN 

PRINCIPAL A 

DAA 

= 

6.599 

IN-LB-S/RAD 

CROSS-COUPLED 

DAB 

= 

-1.312 

IN-LB-S/RAD 

CROSS-COUPLED 

DBZ 

= 

-0. 1294E-03 

IN-LB-S/IN 

CROSS-COUPLED 

DBA 

= 

-1.299 

IN-LB-S/RAD 

PRINCIPAL B 

DBB 

= 

0.3091 

IN-LB-S/RAD 


95TR34-V2 

7.3 Hydrostatic Recess Seal 

The configuration of the hydrostatic recess seal is shown on Figure 50. It is a four-recess seal and 
variable grid was used to get a fine grid definition at the recess extremities where large clearance 
changes occur. Since the entire seal was modeled, a large number of variable grid points must be 
entered. For this example, 137 entries in the circumferential direction were inputted. For 
hydrostatic recesses, the additional input variables are the recess depths and locations and the 
location and sizes of the feeding orifices. 

The land clearance is 0.0005 in. and the recess depth is 0.020 in. One orifice feeds each recess at 
a diameter of 0.010 in. Other important parameters are: 

• Supply pressure = 200 psig 

• Absolute temperature - 1460° R (1000°F) 

• Gas viscosity = 5.35 x 10~ 9 lb-s/in. 2 

• Speed = 0 

• OD pressure = 200 psig 

• ID pressure = 0 psig. 


The output summary is shown on Table 21. Note that the damping is negative indicating an 
unstable seal. Hydrostatic gas lubricated seals or bearings are prone to pneumatic hammer or 
chattering of the opposed surfaces. Generally, recesses should be avoided and gas introduced into 
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the film strictly by the inlet holes. This type of feeding is discussed in the next sample problem. 
The clearance and pressure distributions are shown on Figures 51 and 52, respectively. 


7.4 Inherently Compensated Hydrostatic Seal 

A hydrostatic seal with a series of feed holes that enter directly into the film employs inherent 
compensation which implies that the restrictor area is the curtain area of a cylinder with diameter 
equal to the hole diameter and length equal to the film height at the hole location. The seal 
analyzed had an inside diameter of 2 in. and an outside diameter of 4 in. For this example, holes 
were specified at a radius of 2.4 in and every 10 degrees in the circumferential direction. The hole 
size was 0.010 in. diameter. 

Other parameters are as follows: 

• Clearance = 0.0005 in. 

• Supply pressure = 200 psig 

• Gas constant = 246,900 in. 2 /s 2 /°R 

• Absolute temperature = 1460 °R (1000 °F) 

• Viscosity = 5.35 x 10' 9 lb- s/in. 2 

• Speed = 0 

• OD pressure = 200 psig 

• ID pressure = 0 psig. 


Table 21. GFACE Recess Seal 

-JOURNAL POSITION & FORCES 


MINIMUM CLEARANCE 


= 

0.000500 

IN 

LOAD CAPACITY 


= 

6265. 

LB 

MOMENT ABOUT X-X 


= 

0. 1421E-06 

IN-LB 

MOMENT ABOUT Y*Y 


= 

-0.1740E-06 

IN-LB 

POWER LOSS 


= 

O.OOOOE+OO 

HP 

LEAKAGE AT INNER 


- 

-0.8928E-03 

LB/S 

LEAKAGE AT OUTER 


= 

-0. 2859E-03 

LB/S 

STIFFNESS COEFFICIENTS 




PRINCIPAL 2 

KZZ 

= 

0.3003E+07 

LB/IN 

CROSS-COUPLED 

KZA 

= 

79.56 

LB/RAD 

CROSS-COUPLED 

KZB 

= 

79.73 

LB/RAD 

CROSS-COUPLED 

KAZ 

= 

-0. 1382E-01 

IN-LB/IN 

PRINCIPAL A 

KAA 

= 

0.4779E+07 

IN-LB/RAD 

CROSS-COUPLED 

KAB 

= 

-1.289 

IN-LB/RAD 

CROSS-COUPLED 

KBZ 

= 

-0.3776 

IN-LB/IN 

CROSS-COUPLED 

KBA 

= 

-0.7086 

1 N-LB/RAD 

PRINCIPAL B 

KBB 

= 

0.4779E+07 

IN-LB/RAD 

DAMPING COEFFICIENTS 




PRINCIPAL Z 

DZZ 

= 

-0.1 258E+05 

LB-S/IN 

CROSS-COUPLED 

DZA 

= 

-4.177 

LB-S/RAD 

CROSS-COUPLED 

DZB 

= 

-4.178 

LB-S/RAD 

CROSS-COUPLED 

DAZ 

= 

0.7818E-04 

IN-LB-S/1N 

PRINCIPAL A 

DAA 

= 

-0.1 440E+05 

I N-LB-S/RAD 

CROSS-COUPLEO 

DAB 

= 

0.7198E-02 

IN-LB-S/RAD 

CROSS- COUPLED 

DBZ 

= 

0. 1054E-05 

IN-LB-S/IN 

CROSS-COUPLED 

DBA 

= 

0.3959E-02 

IN-LB-S/RAD 

PRINCIPAL B 

DBB 

= 

-0. 1440E+05 

IN-LB-S/RAD 
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The output summary for this case is shown on Table 22. Notice that the direct damping in the 
axial direction (DZZ) is a positive number and that this configuration is not unstable as is the prior 
problem that incorporated recesses. There is no trapped volume that can cause out-of-phase 
stability problems. 

The clearance and pressure distributions are shown on Figures 53 through 55. 


7.5 Radial Taper Seal 

One of the options in GFACE is the incorporation of a radial taper in the seal face. A tapered seal 
was analyzed. The following parameters were applied: 

• Inner diameter = 1.0 in. 

• Outer diameter =1.5 in. 

• Clearance = 0.0004 in. 

• Pad angle = 360° 

• Specific heat = 1.4 

• Gas constant = 255,488 in. 2 /s 2 /°R 

• Absolute temperature = 520 °R (60 °F) 

• Viscosity = 2.5 x 10' 9 lb-s/in. 2 

• Speed = 30,000 rpm 

• OD pressure = 0 psig 

• ID pressure = 400 psig 

An ID taper was incorporated (see Table 23) that extended half way up the face of the seal, and 
the slope of the taper was 0.1833. Clearance and pressure distributions are indicated on Figures 
56 and 57. 
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Table 22. GFACE Hydrostatic Source Seal 


-JOURNAL POSITION & FORCES 


MINIMUM CLEARANCE 


= 

0.000500 

IN 

LOAD CAPACITY 


= 

6588. 

LB 

MOMENT ABOUT X-X 


= 

0.9286E-09 

IN-LB 

MOMENT ABOUT Y-Y 


= 

-0.8015E-09 

IN-LB 

POWER LOSS 


= 

o 

o 

+ 

UJ 

o 

o 

o 

o 

o 

HP 

LEAKAGE AT INNER 


= 

-0.1 165E-02 

LB/S 

LEAKAGE AT OUTER 


= 

-0. 1727E-03 

LB/S 

STIFFNESS COEFFICIENTS 




PRINCIPAL Z 

KZZ 

= 

0. 1722E+Q7 

LB/IN 

CROSS-COUPLED 

KZA 

= 

410.1 

LB/RAD 

CROSS-COUPLED 

KZB 

- 

410.1 

LB/RAD 

CROSS-COUPLED 

KAZ 

= 

-0.2244E-03 

IN-LB/IN 

PRINCIPAL A 

KAA 

= 

0.2562E+07 

IN-LB/RAD 

CROSS-COUPLED 

KAB 

- 

0.6152E-03 

IN-LB/RAD 

CROSS-COUPLED 

KBZ 

= 

-0.1 138E-03 

IN-LB/IN 

CROSS-COUPLED 

KBA 

= 

0. 1788E-02 

IN-LB/RAD 

PRINCIPAL B 

KBB 

= 

0.2562E+07 

IN-LB/RAD 

DAMPING COEFFICIENTS 




PRINCIPAL Z 

DZZ 

= 

4538. 

LB-S/IN 

CROSS-COUPLED 

DZA 

= 

-0.1395 

LB -S/RAD 

CROSS-COUPLED 

DZB 

= 

-0.1395 

LB-S/RAD 

CROSS-COUPLED 

DAZ 

= 

0.8589E-07 

IN-LB-S/IN 

PRINCIPAL A 

DAA 

= 

0.2322E+05 

IN-LB-S/RAD 

CROSS-COUPLED 

DAB 

= 

-0.2646E-06 

I N-LB-S/RAD 

CROSS-COUPLED 

DBZ 

= 

0.2733E-10 

IN-LB-S/IN 

CROSS-COUPLED 

DBA 

= 

-0.7288E-06 

IN-LB-S/RAD 

PRINCIPAL B 

DBB 

= 

0.2322E+05 

IN-LB-S/RAD 


95TR34-V2 


Table 23. GFACE Inside Radial Tapered Land Seal 


-JOURNAL POSITION & FORCES 


MINIMUM CLEARANCE 


= 

0.000400 

IN 

LOAD CAPACITY 


= 

277.4 

LB 

MOMENT ABOUT X-X 


= 

-0. 1214E-12 

IN-LB 

MOMENT ABOUT Y-Y 


= 

-0.8861E-13 

IN-LB 

POWER LOSS 



0.3356E-02 

HP 

LEAKAGE AT INNER 


— 

0. 1264E-01 

LB/S 

LEAKAGE AT OUTER 


= 

0. 1246E-01 

LB/S 

STIFFNESS COEFFICIENTS 




PRINCIPAL Z 

KZZ 

= 

0.5404E+05 

LB/IN 

CROSS-COUPLED 

KZA 

= 

47.84 

LB/RAD 

CROSS-COUPLED 

KZB 

= 

47.84 

LB/RAD 

CROSS-COUPLED 

KAZ 

= 

-0.3066E-07 

IN-LB/IN 

PRINCIPAL A 

KAA 

= 

0.2056E+05 

IN-LB/RAD 

CROSS- COUPLED 

KAB 

= 

321.0 

IN-LB/RAD 

CROSS-COUPLED 

KBZ 

= 

0.2232E-07 

IN-LB/IN 

CROSS-COUPLED 

KBA 

= 

-321.0 

IN-LB/RAD 

PRINCIPAL B 

KBB 

- 

0.2056E+05 

IN-LB/RAD 

DAMPING COEFFICIENTS 




PRINCIPAL Z 

DZZ 

= 

1.088 

LB-S/IN 

CROSS-COUPLED 

DZA 

= 

-0. 1 101E-03 

LB-S/RAD 

CROSS-COUPLED 

DZB 

= 

-0.1 101E-03 

LB-S/RAD 

CROSS-COUPLED 

DAZ 

= 

0.8441 E- 13 

IN-LB-S/IN 

PRINCIPAL A 

DAA 

= 

0.2047 

IN-LB-S/RAD 

CROSS-COUPLED 

DAB 

= 

-0.1491E-02 

IN-LB-S/RAD 

CROSS-COUPLED 

D8Z 

= 

-0.2549E- 16 

IN-LB-S/IN 

CROSS-COUPLED 

DBA 

= 

0. 1 49 IE- 02 

IN-LB-S/RAD 

PRINCIPAL B 

DBB 

= 

0.2047 

IN-LB-S/RAD 
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GFACE 

FILM THICKNESS DISTRIBUTION 

RAYLEIGH STEP SEAL EXAMPLE 41 


0.D. _ 9.000 IN 

SPEED = 5000-00 RPM 

I -0. = 2-500 IN 


CLEARANCE = 0.000200 IN 
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Figure 44. Sample Problem 1, Rayleigh-Step Seal , Clearance Distribution 
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Figure 45. Sample Problem 1, Rayleigh-Step Seal, Pressure Distribution 
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PRESSURE DISTRIBUTION 

SPEED = 5000.00 RPM 
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GFflCE 


FILM THICKNESS DISTRIBUTION 

SINGLE TRPERED 

LRND SERL 


O.D. - 

4.505 IN 

SPEED = 50000.00 RPM 

I -D ■ = 

3-793 IN 


CLERRRNCE _ 0 

000200 IN 
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Figure 48. Sample Problem 2 y Tapered Land Seal , Clearance Distribution 
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Figure 49. Sample Problem 2, Tapered Land Seal, Pressure Distribution 
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GFRCE FILM THICKNESS DISTRIBUTION 

RECESS SERE , 4 PRDS . JOINED BOUNDARIES 

0-D. - 8.000 IN SPEED = 0-00 RPM 

I -D- = 4 -000 IN 

CLERRRNCE _ 0-000500 IN 



•^o 






Figure 51. Sample Problem 4 , Four-Recess Seal, Clearance Distribution 
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Figure 52. Sample Problem 4, Four-Recess Seal, Pressure Distribution 
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GFRCE 


FILM THICKNESS 

DISTRIBUTION 

HYDRGSTRTIC 

SOURCE SERL 




0.D. = 
I -D - - 
CLERRRNCE - 

8.000 IN 
4-000 IN 
0.000500 IN 

SPEED r 

0.00 RPM 

■ 




_ .500E-03 




\ 

- 




-^^360 
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Figure 53. Sample Problem 6, Inherently Compensated Seal , Clearance Distribution 
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Figure 54. Sample Problem 6, Inherently Compensated Seal , Pressure Distribution 
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OFFICE 

HYDRGSTRTIC SOURCE SEAL 
0.D. ^ 8.000 IN 

I-D- = 4.000 IN 

CLEARANCE = 0-000500 IN 


PRESSURE DISTRIBUTION 

SPEED = 0.00 RPM 
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Figure 56. Sample Problem 7 , Radial Taper Seal, Clearance Distribution 
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8.0 VERIFICATION OF CODE GFACE 


8.1 Tilted Slider 

Etsion and Fleming 191 found that maximum load capacity of a slider occurs when it is tilted about 
the trailing edge and they produced analytical data for varying compressibility numbers, sector 
angles and radius ratios. Table 24 below shows comparisons for a 45 degree sector. 


Table 24. Comparisons of GFACE with Etsion and Fleming Load Capacity 


A 

s 

c 

degrees 

N 

rpm 

W 

GFACE 

W' Etsion 
& 

Fleming 

A% 

1 

2.5 

.0143 

133.69 


0.0047 

1.27 


2.5 

.0143 

1,336.9 


0.0460 

1.30 

25 

3.0 


3,342.25 

0.1086 

0.1080 

0.56 


3.5 


6,684.5 

0.1897 

0.1850 

2.5 


4.5 

.0258 

13,369 

0.3014 

0.2950 

2.17 


where: 


A = 


6(Ltcor l 


£ = 


C r 2 


_ b x 2 . 


W' = 


2W 


P ~ 




lb — s 

|i = absolute viscosity , = 1.75 x 10“ 9 — 

in. 

, _ , rad 

co = shatt speed, 

s 

N = shaftspeed, rpm 
W = load, lb 

P a = ambient pressure, = 14.7psia 
r 2 = outside radius, = 2.0 in. 

Tj - inside radius, = 0.6 in. 

Q = tilt angle, radians 
£ = tilt angle, degrees 

h 2 = trailing edge film thickness = 0.0002 in. 

The comparative values of non-dimensional load, W are within 3%. Comparisons were also 
made for the Power Loss Coefficient (PLC). Results are indicated on Table 25. 
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Table 25. Comparison of GFACE with Etsion and Fleming Power Loss 


A 

£ 

PLC 

GFACE 

PLC 

Etsion 

1 

2.5 

49.4 

12 

10 

2.5 

12.86 

12 

25 

3.0 

12.987 

13 

50 

3.5 

14.193 

15 

100 

4.5 

16.15 

17 


PLC = where: 

Wcoh 2 

F = Power Loss, in. -lb / s 

The correlation is excellent except for A =1. The friction levels are on a precipitous upward 
curve at this value of A and it is very easy to be off by a substantial amount. The curve fitting of 
the Etsion data could easily produce significant error. It is the author’s inclination that the 
GFACE values are accurate. 


8.2 Ausman Rayleigh-Step Analysis 

In 1961 Ausman produced theoretical results for a Rayleigh step sector [10] for varying radius 
ratios, number of pads and compressibility parameter A. Comparative results with GFACE are 
indicated on Table 26. 


Table 26. GFACE Comparison with Ausman Rayleigh-Step Pad 


A 

n 

W 

W 

GFACE 

W’ 

Ausman 

A % 

10 

8 

0.7539 

0.0456 

0.046 

-.87 

20 

7 

1.821 

0.0957 

0.103 

-7 

40 

7 

3.701 

0.1946 

0.219 

-11 

80 

6 

7.766 

0.3479 

0.397 

-12 

160 

6 

11.71 

0.5246 

0.572 

-8 


6 poor 2 

php 


n ~ number of pads; 


r i 


= 0.5 


r 2 


Variables remain as defined above. The correlation as indicated is fair, but not as good as 
expected. The Ausman analysis was done in 1961, when contemporary numerical analysis was 
unavailable. Therefore, it is believed that the GFACE data is more accurate. 
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Further cases were compared with information in Gross 1111 for a step pad configuration with the 
following geometry and operating conditions: 


Shaft speed, N 
Ambient pressure, P a 
Viscosity, p, 

Radius ratio, R 
Outside radius, 

Inside radius, r j 
Pad Angle, (3 
Step Angle, (3] 

Land film thickness, h 2 


36,000 rpm 
42.6 psia 
5 x 10~9 lb-s/in.^ 
0.6 

1.75 in. 

1.05 in. 

38° 

150 

0.0005 in. 


Comparative results are shown on Figures 58 and 59. The comparisons are reasonable and 
discrepancies are probably due to the inaccuracies of the numerical methods utilized by Gross that 
appears to overestimate load capacity and underestimate power loss. Thus, the application of 
GFACE would provide conservative designs as compared to Gross. 


8.3 Internal Verification 


There are a number of internal checks that can be made to assure numerical correctness and 
accuracy. The first case is a recess hydrostatic bearing without rotation. Then the mass flow 
entering the recess from an external source should equal the mass flow exiting the pad. A 35 
degree pad was examined with variable grid geometry and a recess specified. Four flow paths 
were specified at the perimeters of the pad. The computer output indicated the following: 


Flow exiting from pad: 

• Flow path 1 = 

• Flow path 1 - 

• Flow path 1 = 

• Flow path 1 = 

• Total = 


-0.1968 x 10 " 4 lb/s 
0.3519 x 10 - 4 lb/s 
-0.6005 x 10 - 4 lb/s 
0.8530 x 10 - 4 lb/s 
2.0022 x 10 - 4 lb/s 


The negative number implies that the flow is in a direction opposite to that of the positive sign 
convention. In this case all flows are exiting the pad and the total outflow is the absolute sum. 


NAS A/CR— 2004-2 1 3 1 99/VOL2 


115 



Load (lbs) 


Step Validation, Gross 


40 

35 

30 

25 

20 


10 

5 

0 

0123456789 10 

Step Height x 10**4 

95TR34-V2 

Figure 58. Validation , Step Pad, Load vs. Step Height 
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Step Validation, Gross 



Step Height x 10**4 

95TR34-V2 


Figure 59. Validation , Step-Pad , Load/Power vs. Step Height 
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The inflow is the total flow going through the orifice feeding the recess. The following equation 
applies: 


= 386.4 A 0 C d G,P s 


( v 3 


v p »y 




Y-i 

Y 


v p s y 



7td, 


nO.OlO' 


A„ - orifice area = 

0 4 4 

C D = Coefficient of Discharge = 0.9 

P r _ recess pressure _ 131 + 14.7 

P s supply pressure 200 + 14.7 


= 7.854 x 10 


-5 


= 0.678621 


2y 


° l V g c 0(y-i) 

where y = ratio of specific heats = 1 .4 


in. 


G = gas constant =246,900 . 
c s 2 -° R 

0 = absolute temperature = 1460° R 
f = 386.4x7.854xl0” 5 x0.9xl.3935xl0 _4 x214.7 


x< (0.678621)14 


0.4 


1- (0.678621) w 


IK 

= 2.006x10"* — 

S 


The correlation between the exit and inlet flows is excellent. 

A check on the stiffness can be made, at zero excitation frequency, by checking the stiffness 
values against incremental displacements in the three degrees of freedom. A Rayleigh step pad 
was examined. Pertinent parameters are as follows: 

Outside diameter = 4 in. 

Inside diameter = 2.5 in. 

Clearance = 0.0002 in. 

Step height = 0.004 in. 

Specific heat =1.4 

Gas constant = 246,900 in 2 /s 2 /°R 

Absolute temperature = 1460 °R (1000 °F) 

Operating Speed = 5000 rpm 
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Perimeter pressure = 0 psig 
Pad angle = 35° 

Step angle = 25° 

Step width = 0.43 in. (centrally located). 

Table 27 shows the results of the displacements on Forces and Moments. 

Table 27. Forces and Moments as a Function of Displacements 
(3z = 0.0001 in., 3a - 0.0001 degrees., 3ft - 0.0001 degrees) 


Displacement 

0 

Z 

a 


Axial Force, F 7 lb 

4.861 

4.531 

4.854 

4.869 

Moment about x- 
axis, M^y in. -lb 

3.838 

3.571 

3.832 

3.844 

Moment about Y- 
axis, M w in. -lb 

-6.946 

-6.479 

-6.936 

6.957 


The stiffness is computed by subtracting the equilibrium values from the displaced values and 
dividing by the displacement. Since a positive displacement results in a negative force differential, 
the quantities are multiplied by -1 to maintain the accepted stiffness sign convention. Table 28 
below compares the computed values of stiffness as calculated by the code against those manually 
computed by the above procedure. The top value is computed by exercising the stiffness option 
of the code and the bottom value is the manually computed number. The corroboration is 
excellent. Additional runs were made with full 360° pads with similar results. 


Table 28. Stiffness Comparisons Between Automatic and Manual 
Computed Values; Rayleigh-step Pad 
(3z - 0.0001 in., 3a = 0.0001 degrees., 3ft = 0.0001 degrees) 


K zz 

K az 


34,570 

27,930 

-48,980 

33,000 

26,700 

-46,700 




K za 

K r/(7 


39,900 

31,750 

-57,150 

40,107 

34,377 

-57296 




KHHHi 



-44,200 

-35,840 

62,710 

-45,836 

-34,337 

63,025 
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